NAVAL POSTGRADUATE SCHOOL 
IVIonterey , California 




THESIS 



EFFECTS OF FADING AND DATA MODULATION 
ON NONCOHERENT M-SEQUENCE ACQUISITION 

SCHEMES 

by 

Patrick J. Vincent 
March 1992 

* 

Thesis Advisor Professor Alex W. Lam 



Approved for public release; distribution is unlimited. 



T254592 



Unclassified 



security classification of this page 



REPORT DOCUMENTATION PACE 


la Report Security Classification Unclassified 


lb Restrictive Markings 


2a Security Classification Authority 


3 Distribution Availability of Report 

Approved for public release; distribution is unlimited. 


2b Declassification Downgrading Schedule 


4 Performing Organization Report Number(s) 


5 Monitoring Organization Report xNumber(s) 


6a Name of Performing Organization 
Naval Postgraduate School 


6b Office Symbol 
(if applicable) 32 


7a Name of Monitoring Organization 
Naval Postgraduate School 


6c Address (dry, state, and ZIP code) 

Monterey, CA 93943-5000 


7b Address (city, state, and ZIP code) 

Monterey, CA 93943-5000 


8a Name of Funding Sponsoring Organization 


8b Office Symbol 
( if applicable ) 


9 Procurement Instrument Identification Number 


8c Address (city, state, and ZIP code) 


10 Source of Funding .Numbers 


Program Element No Project .No Task .No Work Unit Accession .No 


11 Title (Include securirv classification) ETFECTS OF FADING AN D DATA MODULATION ON NONCOFIFRFNT 

m-sfouence acouisition schemes 


12 Personal Author(s) Patrick J. Vincent 


13a Type of Report 13b Time Co\ered 

Master's Thesis From To 


14 Date of Report (year, month, day) 15 Page Count 

March 1992 110 



16 Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official policy or po- 
sition of the Department of Defense or the U.S. Government. 



17 Cosaii Codes 


Field 


Group 


Subgroup 















1$ Subject terms {continue on reverse if necessary and identify by block number) 

sequential acquisition, spread spectrum, communications 



19 Abstract (continue on reverse if necessary and identify by block number) 

In direct -sequence spread -spectrum systems, successful communications require phase synchronization of the m-sequence 
in the incoming signal with a locally generated m-sequence at the receiver. Many acquisition schemes which extract the phase 
of an incoming m-sequence have been studied, but most of them assume coherent demodulation (which is usually not 
available during acquisition) and 'or independent samples (which introduce a loss in the effective signal to noise ratio (SNR)). 

This thesis investigates the performance of two m-sequence acquisition schemes in the presence of fading and data mod- 
ulation. A fixed sample size test and a truncated sequential test are studied without the usual assumptions of coherent de- 
modulation or independent samples. The effects of fading and data modulation on our schemes' probability of false alarm, 
probability of detection and test length are thoroughly explored. We find that channel fading in effect induces a loss of signal 
SNR, but the desired power of the tests can be restored by suitable adjustments in the decision processor. We find that the 
effects of data modulation are less severe, but more problematic to correct. 



20 Distribution Availability of Abstract 

(3 unclassified unlimited D same as report D DTIC users 



21 Abstract Security Classification 
Unclassified 



22a Name of Responsible Individual 
Professor Alex W. l.am 



22b Telephone (include Area code) 

(408) 646-3044 



22c Office Symbol 

EC/LA 



DD FORM 1473,84 MAR 



83 APR edition may be used until exhausted 
All other editions are obsolete 



security classification of this page 

Unclassified 



1 



Approved for public release; distribution is unlimited. 



Effects of Fading and Data Modulation 
on Noncoherent M-Sequence Acquisition Schemes 

by 

Patrick J. ^uncent 
Lieutenant, United States Na\y 
B.S., Polytechnic Institute of New York, 1984 

Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 

from the 

NAVAL POSTGIUADUATE SCHOOL 
March 1992 



ABSTRACT 



In direct-sequence spread-spectrum systems, successful communications require 
phase synchronization of the m-sequence in the incoming signal with a locally generated 
m-sequence at the receiver. Many acquisition schemes which extract the phase of an 
incoming m-sequence have been studied, but most of them assume coherent demodu- 
lation (which is usually not available during acquisition) and'or independent samples 
(w’hich introduce a loss in the effective signal to noise ratio (S\R)). 

This thesis investigates the performance of two m-sequence acquisition schemes in 
the presence of fading and data modulation. A fixed sample size test and a truncated 
sequential test are studied without the usual assumptions of coherent demodulation or 
independent samples. The effects of fading and data modulation on our schemes' 
probability of false alarm, probability of detection and test length are thoroughly ex- 
plored. We find that channel fading in effect induces a loss of signal SNR, but the de- 
sired power of the tests can be restored by suitable adjustments in the decision processor. 
We find that the effects of data modulation are less severe, but more problematic to 
correct. 
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I. INTRODUCTION 



This thesis explores the elTects of fading and data modulation on sequential 
pseudonoise (PN) acquisition schemes. We begin by providing a brief introductor>‘ 
overview of spread spectrum communication systems. The purpose of this introductory 
chapter is twofold. The inclusion of this material enhances the degree to which this 
thesis is self-contained. Additionally, the notation introduced in this chapter will be used 
in subsequent sections. 

A complete presentation of the theorj' and applications of spread spectrum com- 
munication systems would fill several volumes. The intent here is neither to be com- 
prehensive in scope nor specific in detail. We present, in this chapter, only a thumbnail 
sketch of the subject, including only those details necessar>’ to establish a foundation for 
the following chapters. 

Spread spectrum signals allow us to [Ref 1 : pp. 800-801): 

1. overcome or minimize the elTects of interference, whether it be intentional (e.g., 
jamming), unintentional (e.g., interference from other channel users) or self-induced 
(e.g., multipath). 

2. hide a signal by "burying" it below the background noise level, thus concealing the 
transmission from unintended listeners and eavesdroppers. 

3. achieve message privacy. 

In a conventional communication system employing amplitude modulation (A.M) 
or phase modulation (PM), the communication signal bandwidth is approximately equal 
to twice the rate of data transmission. In spread spectrum (SS) systems the bandwidth 
of the communication signal is much greater than the data rate. In all SS systems, ran- 
domness is used to make interception and copying more formidable. 

Prior to 1980, the field of spread spectrum communications fell almost entirely 
within the military domain. Techniques were developed to allow communications which 
were not detectable by adversaries and to combat the disruptive effects of jamming by 
an opponent. In 1983 the Federal Communication Commission (FCC) opened, for the 
first time, three frequency bands for commercial spread spectrum use [Ref 2J. Cellular 
and indoor commercial networks using spread spectrum techniques are now becoming 
a reality. 

In the following w'e will discuss several typical spread spectrum techniques. 
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A. TIME HOPPING (TH) 

In a time hopping scheme, we transmit the communication signal at randomly cho- 
sen times for random burst durations. If ri is the fraction of time the transmitter is active 
(7 < 1 ) and r, is the steady data rate from the user, then , the transmission rate when 
the transmitter is activated, is given by 







( 1 . 1 ) 



Note that the bandwidth of the signal that gets transmitted is proportional to r* (which 
may be much larger than ) so we have distributed, or spread, the signal over a wider 
frequency band. 

The primar>' advantage of TH is an improvement in the Signal to Jamming Ratio 
(SJR). Let S and J be the average signal power and average jamming power, respec- 
tively, at the receiver input. Then the SJR for a non-SS system is given by 

SdR = j- ( 1 . 2 ) 

For a TH system, let S be the input power at the receiver when a burst transmission 
occurs. Then S = St] and the SJR during burst reception is 



A 




Thus we see that the SJR in a TH system is better than that for a non-SS system. The 
amplification factor (= -^) is called the processing gain (PG). 

The improvement in SJR in a time hopping system is countered by several disad- 
vantages: 

1 . a data buffer must be employed in both the transmitter and receiver 

2. the TH transmitter must radiate at higher power levels 

3. the TH receiver must resynchronize on each new burst 

B. FREQUENCY HOPPING (FH) 

As previously noted, in all SS systems randomness is used to make the communi- 
cations system more immune to signal detection and jamming. In a system employing 
frequency hopping, the transmitter carrier frequency changes value (hops) as a function 
of time. If wc let be the amount of time the carrier remains at each of the hopping 
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frequencies, then the bandwidth of the carrier is always constant and equal to 2//* = b 
Hz. If B is the total range of frequencies we are assigned to hop over, then the number 
of distinct carriers used is equal to Bib. 

A FH receiver despreads the received signal by multiplying it by a locally generated 
hopping signal. Note that the hop sequence of the receiver's local oscillator must match 
the hop sequence of the transmitter in order to recover the signal. 

As with time hopping, a major advantage of FH systems is an improved immunity 
to jamming. In a FH scheme, the jamming power is now spread over a bandwidth of 
B Hz. The signal bandwidth is approximately equal to b Hz. Although the entire signal 
power is present after dehopping, the amount of jamming power present after dehopping 
is JfbiB). Thus, in a FH system. 




The processing gain is thus given by B'b. 

In many military' applications we would like to deny an adversary the knowledge 
that we are transmitting a signal. That is to say, besides making it difficult for an ad- 
versar}' to "break" our code once he intercepts our signal, we would like to make it dif- 
ficult for him to detect that there exists a transmitted signal at all. (Obviously a signal 
that is not detected can not be exploited.) To achieve a low probability of detection in 
a communications scheme, the signal power level must be reduced below the noise power 
level. In other words, we need to "bur\" the signal in the noise (or hide it behind another 
more powerful already-present signal) so that the adversary' will "see" only noise. Such 
a scheme is termed a Low Probability of Detection (LPD) scheme. 

A major disadvantage of the frequency hopping approach lies in the fact that the 
signal can not be reduced below the noise level; in fact, the signal must be maintained 
significantly above the noise level. As a result, frequency hopping schemes are used to 
provide anti-jamming, not low probability of signal detection. 

In some communication environments, multiple transmission paths may arise be- 
tween the transmitter and receiver. This phenomenon, termed multipath, often degrades 
communications since the signals arriving at the receiver via secondary paths 
destructively interfere with the signal on the primary' link. In a FH receiver only the 
energy at a particular frequency is passed through the first stage filter at any interval of 
time. If the multipath signals arrive at the receiver with a sulTicient time delay, then the 
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receiver may have "gone on lo the next hop" by that time. In such instances the multi- 
path signals are effectively ignored. This is called the multipaih immunity of frequency 
hopping. 

Finally, the utility of the frequency hopping method is limited by the speed with 
which the electrical components can alter frequency without generating excessive noise. 
As of 1990, frequency oscillators were able to var}’ frequencies at rates of 1 million hops 
per second over a 20 .MHz bandwidth [Ref 3). 

C. DIRECT SPREADING (DS) 

1. General 

In a direct spreading scheme, the bandwidth spreading is accomplished on a 
bit-by-bit basis. We multiply the data signal Vi(/) by a wideband voltage signal vft). 
Since multiplication in the time domain corresponds to convolution in the frequency 
domain, the resulting signal's bandwidth is spread. We can then use a matched filter or 
correlator in the receiver to detect the transmitted signal. 

If a matched filter is used in the receiver, an e.xpression can be obtained for the 
system's processing gain. Let E represent the energy per bit in the input signal to the 
receiver, let T represent the bit duration, let IF be the bandwidth of the spread signal and 
let NJ2 be the two-sided power spectral density of the noise. Then the input signal 
power is given by 




(1.5) 



The input noise power is given by 




( 1 . 6 ) 



and the input SNR is 




(1.7) 



A result from matched filter theory states that the output signal to noise ratio is given 



Combining (1.7) and (1.8) yields 



S.V/?, = SA'/?J2riF) (1.9) 

Thus 2TIV is the processing gain for direct spreading using a matched filter. To in- 
crease the processing gain we usually increase IV. Note that there is no theoretical limit 
to the degree of improvement achievable. 

It is a classic myth to believe that direct spreading improves the output SNR in 
receivers where thermal noise limits operation. This is not the case. If the receiver's 
performance is limited by thermal noise, then employing spread spectrum will not help. 
Looking at (1.9), we know that utilizing spread spectrum increases JV and hence, we may 
be tempted to conclude, the output signal to noise ratio, SjVR^ , must increase. How- 
ever, increasing li' increases the noise bandwidth of the receiver, and hence increases the 
input noise power since 



N-=kTW (l.IO) 

where k is Boltzmann's constant and T is the effective noise temperature. The increase 
in A’, causes a decrease in The decrease in S.\'R,„ counterbalances the increase in 

the (2TIV) term. 

In jamming limited receivers A' = 7 =constant. In this case we do get improve- 
ment by employing direct spreading. 

Direct spreading alTords us the ability to manage the power in the transmitted 
signal. We can distribute the power over any frequency interval without penalty. We 
can distribute the power over a wide frequency band so that the signal power density 
becomes less than the noise power density. In this way we can "bury" the signal in the 
noise and achieve covert communications. 

To summarize, in the direct spreading technique, we represent each bit of dura- 
tion T with a wideband signal v,(0 duration T and bandwidth W. 

2 . The Wideband Signal 

The value of the direct spreading scheme depends critically on the choice of 
v,(r), the wideband signal. The signal must have a good autocorrelation function (ACF) 
since we want a good matched filter output. A "good" ACF is one with a narrow , high 
magnitude main lobe and small side lobes. Such an ACF allows the use of a follow-on 
threshold device with minimum decision error. Additionally, v,(/) should have a 
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buildable matched filter and should have a near-white (flat) spectrum to allow for covert 
communications. 

3. Analog Direct Spreading 

In analog direct spreading, the wideband signal v,(/) is analog. The beauty of 
analog direct spreading is that the communication signal can be directly transmitted; 
there is no need for a separate modulating carrier. Perhaps the most commonly used 
analog spreading signal is the CHIRP signal. In this case v,(t) varies linearly in fre- 
quency from some initial frequency,/, to some final frequency,/, over the interval T. 
The formula for the spreading signal is 

Vjf(/) = A cos 

To send the bit "1" we may start at/ and terminate at/ (called UP-CHIRP if/>/). 
To send the bit "0“ we start at / and terminate at / (called DOWX-CHIRP if/>/). 
Using Fresnel Integrals it can be shown that the bandwidth of the spread signal is ap- 
proximately equal to / — /. Additionally, the spectrum is appro.ximately constant in the 
interval (/J,/) and zero outside this interval. 

The ACF of the analog CHIRP signal is quite good. Figure 1 on page 7 shows 
one period of the autocorrelation function of v.(t) given by (1.11) for the case of 
2TW = 100, r= Isec and / = 50 Hz. The matched filter for CHIRPED signals is the 
surface acoustic wave device (SAW) (Ref 4]. 

4. Digital Spread Spectrum 

In digital direct spreading, the wideband signal v/t) is digital. Specifically, for 
v,(/) we choose a binary sequence. Instead of transmitting a data bit (say a "1") we 
transmit a binary sequence (say, "1100011"). Each bit in the binarj' sequence that 
comprises v^(/) is called a chip. As of 1990, the maximum achievable chip rate was 50 
million chips per second [Ref 3). 

The digital spreading signal is created by using the output of an «-stage shift 
register with feedback connections. An «-stage shift register can take on, at most, 
2" — 1 different states (depending on the feedback connections), so the maximum length 
sequence attainable from the shift register is also equal to 2" —1. These maximum length 
sequences, called m-sequences, have the best ACF attainable for sequences of their 
length. The theory’ of feedback shift register sequences is a well-developed subject, and 
the interested reader is referred to (Ref 5 : Chapter 7]. 
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Figure 1. Autocorrelation Function of CHIRP 

The ACFs of all m-scqucnces have the same form. The peak value is propor- 
tional to 2" —1. The sidelobe value is proportional to -1. A normalized m-sequencc ACF 
is shown in Figure 2 on page 8. 

Sending data using a spreading sequence (e.g., an m-sequence) is called Direct 
Sequence Spread Spectrum (DSSS). To send a "1" we can send one period of the m- 
sequence. To send a "-1" w'e can send one period of the inverted m-sequence. For the 
remainder of this thesis we limit our consideration to DSSS systems. 

5. Matched Filters vs. Correlators 

Matched filters for binarj' sequences usually consist of some variant of a tapped 
delay line (realized as a tapped shift register). The received signal is clocked into a shift 
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Figure 2. Autocorrelation Function of an M-Sequence 

Source: Couch, L.W. II, Digital and Analog Communication Systems , Macmillan, 1990 

register (of length equal to one m-sequence period) in the receiver. The contents of each 
individual register are tapped, inverted as appropriate and summed to form a discrete 
approximation to the ACF. 

Unfortunately digital matched filters have many deficiencies. A large number 
of taps are required for desired values of processing gain, making the hardware imple- 
mentation (on an integrated circuit (IC) chip) difficult if not impossible. Additionally, 
the chips arc "analog" to the extent that they arc buried in noise, and so digital hardware 
does not readily lend itself However we can still use the ACF property of m-sequences 
and replace the matched filter with a three-port correlator. Unfortunately, but 
unavoidably, data recover)’ with a correlator requires a receiver-based locally generated 
matching m-sequence, phase synchronized with the received signal's m-sequence. Of 
course generating the local m-sequence is of no difficulty, but efficiently phase-locking 
it with the received signal has been, and continues to be, a major research concern. 
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II. NONCOHERENT SEQUENTIAL ACQUISITION OF M-SEQUENCES 

In this chapter we consider schemes used to acquire m-sequences in direct sequence 
spread spectrum (DSSS) communication systems. l 

A. INTRODUCTION 

In order to recover the data from a DSSS communications signal, the signal must 
be despread, i.e., the data must be extracted by removing the m-sequence. Because 
matched filters are impractical in this application, the autocorrelation property of m- 
sequences must be exploited. One way to remove the m-sequence (and preserve the 
data) is to multiply the incoming signal by a locally generated identical m-sequence 
which is phase synchronized with the incoming signal's m-sequence. Thus, a primary 
function of the receiver is determination of the phase of the m-sequence in the incoming 
signal. The phase is determined by correlating the incoming signal with the m-sequence 
generated in the receiver. The local sequence's phase is adjusted until synchronization 
is achieved. Aligning the two m-sequences is accomplished in two stages. The first 
stage, called acquisition, brings the m-sequences into coarse alignment. The following 
stage, called tracking, brings the sequences into precise alignment. 

In this thesis we consider a noncoherent acquisition scheme. {Noncoherent, in this 
case, refers to the modulating carrier, not the m-sequence.) Although coherent schemes 
are simpler and more widely studied, they are not practical because m-sequence acqui- 
sition is usually performed prior to recovering the modulating carrier. 

To acquire the incoming signal, we must examine, using some methodology, the 
various phases that the m-sequence can assume. One methodology would be to simul- 
taneously examine all possible m-sequence positions by utilizing, in the receiver, a sep- 
arate correlator for each potential phase. Such parallel schemes are characterized by 
short acquisition times but complex hardware requirements. On the other hand, we can 
examine each possible m-sequence phase one by one to determine if the local sequence 
is aligned with the incoming sequence. Such serial schemes are characterized by long 

1 This chapter, which presents the model and framework to be used in the subsequent chap- 
ters, is adapted from |Ref. 6 ), a paper that served as the origina] motivation for this research, "niis 
chapter, at various times, uses the words, assumptions, notations, equations, methods of organiza- 
tion and conclusions originally propounded in (Ref. 6). This material is adapted and used with the 
permission of the authors. 
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acquisition times but simplified hardware requirements. In this thesis we consider a se- 
rial acquisition scheme. 

B. NONCOHERENT ACQUISITION SCHEMES 

A block diagram of the receiver's acquisition system is shown in Figure 3 on page 
11. The receiver's input signal is comprised of the transmitted signal and corrupting 
additive white Gaussian noise (AWGN) with two-sided power spectral density NJ2. The 
input signal at the receiver, assuming no modulation or doppler shifting, is 



where is the signal amplitude, a(t) is the m-sequence signal waveform with phase 
/AT, (/ taken to be an integer without loss of generality), T, is the chip interval, A is the 
amount that the phase of the local m-sequence is altered during the acquisition process, 
and B are the frequency and phase of the carrier, and n(t) is the AWGN. The locally 
generated m-sequence waveform is described as a[i + (j+ y)ATJ where j is an integer and 
I >’ I < 0.5. The incoming communication signal is multiplied by the locally generated 
signal and then noncoherently demodulated. The result Y„ is checked by the decision 
processor to determine if the local and incoming m-sequences are lined up to within 
Ar,/2. If they are not aligned, the local m-sequence phase is altered by AT, and the 
process repeats. If they are aligned, then j must equal / and the tracking circuit is then 
employed to reduce y to zero. 

Using the notation in Figure 3, and neglecting double frequency terms, we have 



r{t) = Aga{l -t- /ATJ cos((Dp/ -f 0) -b «(/) 



( 2 . 1 ) 




0 




0 



where 




0 



(2.4) 
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Figure 3. Noncoherent Serial Acquisition System 

= f + (/■ + y)Ar^) sin(cOof)i/i (2.5) 

•’o 

are independent, zero mean Gaussian random variables with variance a\ = nT,NJA and 
where 



s„ = 4- f + V+ y^dt (2.6) 

‘c Jq 

For a fixed n, X,^ and the inphase and quadrature components, are independent 
Gaussian random variables, each with variance al and with means 
{AJ2)T^S„co$ 6 and {AJ2)T^S„ sm 6 respectively. Note that X,j, and X^„ (and similarly 
X,^ and X^_„) are not independent for n ^ m,. The test statistic for determining alignment 
is 



n - 4 + P’) 

The random variable Y„ has a probability density function that is noncentral Chi-squared 
with two degrees of freedom [Ref 1 : pp. 25-29], given by 
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( 2 . 8 ) 
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where cos + [(.4„,'2)r,S„ sin 0 ]^ = (/if/4) 7 ^S„^ and where /.(.) is the 

modified Bessel function of order zero. 

Using Y„, the decision statistic, the acquisition scheme determines if the two m- 
sequence phases are aligned within ATJ2, i.e.,j=i=0 (we let / =0 to simplify the nota- 
tion), or if they differ by one or more chips, i.e., | (/ + y) — / 1 AT^ > 7"^ - In other words, 
the decision processor decides between the following hypotheses: 

(non-alignment) : |y -I- y | > -^ and | y | < Yi 



//, (alignment) : 7 = 1 = 0 and | y | < Yi 



(2.9) 



Note that the m-sequence phases may be such that our situation falls between //„ and 
//,. In this case Yi < I 7 + y | < 

The parameter takes on different values under //. and //,. We designate and 

as the worst-case values of under H, and //, respectively. Note that > 0. 

Using these worst-case parameter values, the likelihood ratio for our test becomes 
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We now describe the three decision schemes that we will use to resolve between the 
two hypotheses. 

1. Fixed-Sample-Size (FSS) Scheme 

In a fixed-sample-size (FSS) decision scheme, the length of integration is fixed 
a priori and a decision is made based on the resulting test statistic, F„. If the integration 
length is from r = 0 to r = MT^, the FSS test can be described by 



FSS: 



■^mCv.w) 



> T say Hi 
< T say //q 



( 2 . 11 ) 



where t is the threshold. Since the likelihood ratio is a monotonically increasing func- 
tion of the variable (Ref. 6 ], the FSS test is equivalent to the following test: 
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= AX/(t) say H, 
^■^'= AX/W say Hq 



(2.12) 



where A;'(.) is the inverse function of A„(.). 

2. Sequential Probability Ratio Test (SPRT) 

The sequential probability ratio test (SPRT) consists of testing the likelihood 
ratio against two thresholds for n = 1,2,3-" until one of the thresholds is exceeded. The 
SPRT is optimal for independent Guassian statistics. The length of integration increases 
by one chip each time n increases by one. The threshold that is reached first (upper or 
lower) determines the hypothesis that is accepted (//, or Ha). The test is described by 

j. > A say //, 

SPRIi A„0g \ <B say//o (2-13) 

^ otherwise, continue to next n 

Since Y„ and are not independent random variables we can not write A„(y„) as a 
product of independent random variables. However, we can place bounds on the Type 
I error (false alarm probability) and the Type II error (miss probability) by using Wald's 
inequalities [Ref 7]: 



B B 

^ — , and B>- (2.14) 

where a and fi are the resulting false-alarm and miss probabilities, respectively. If the 
average test length is large, the inequalities can be approximated by equalities. To make 
the real-time SPRT implementation more practical, we can rewrite the test as 

. >A{n) =a;\A) say//, 

SERI: yv ] <B(n)=A;\B) say //q (2.15) 

^ otherwise, continue 

where the thresholds, which are functions of «, can be pre-computed. 

3. Truncated Sequential Probability Ratio Test (TSPRT) 

When using the SPRT it is possible that, for certain m-sequence phase dispari- 
ties, the decision scheme may take an excessively long time to decide between //j and 
//,. Such an incident is especially likely when the difference between the phases of the 
two m-sequences is such that it corresponds to a state between H, and Hg. To avoid a 
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very long test we impose an upper bound on the test length by truncating the SPRT 
(hence the name TSPRT) at a predetermined test duration, and conducting an FSS test 
at that point. The test is described by 







A 

^ > A say //, 




TSPRT: 


if A7 ^ rt. 


A;,(y«) j <B say//o 








^ otherwise, continue 






if/7 = «. 


r >T sav//, 

A„0g = A;te) ( . ■ 

[ <T say//o 


(2.16) 



The test is truncated at n = n (provided it has not already terminated). Again we use the 
monotonicity of A„(.) and rewrite the test as 

. >A{n) = A~\A) say//i 
TSPRT: ifn<n, < <S(«) = A~'(5) say Hq 

^ otherwise, continue 



ifrt = «, w 



>A;’(?) say//, 
<AJ'(t) say Hq 



(2.17) 



A key parameter in determining the utility of an acquisition scheme is the time 
that the test takes to acquire the incoming m-sequence. Generally, for a given false- 
alarm and miss probability, the test with the smallest sample size («) will reach acquisi- 
tion the soonest. The initial phase disparity between the incoming and local 
m-sequences can be regarded as a uniform random variable distributed over the sequence 
period. The hypothesis is, therefore, the most probable and so the expected test 
length under //„ contributes most significantly to the value of the test acquisition time. 
Hence, it is best to minimize the value of the expected sample size under while keep- 
ing the expected sample size to within reasonable values when does not apply. 

C. DESIGN OF DECISION PARAMETERS 

In this section we determine the values of the thresholds to be used in the three ac- 
quisition schemes. 

Denote the chips of the m-sequence by c*, ^ = ...,0,1,2,..., where c* = ±1, and let 
the period be A' = 2" — 1, so that c**;^ = c* for all k. The m-sequence spreading 
waveform can then be written as 
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(2.18) 



^ c/j^r^i - kT,) 

k=—oo 

where pT^{t) is a rectangular pulse of width T, and height 1, i.e., prjii) = 1 for 0 ^ / < 7, 
and it is zero elsewhere. Letting /=0, we obtain expressions for S„ by combining (6) and 
( 18 ): 



S 



n 



n-\ 



f - I y I A)Cjfc + ( I y I A)C;^+5gn(y)] 

I k=0 

1 n-\ 

^C^[(l - + 6^^ +J+;] 



k=0 



under //, 
under Hq 



n-\ 

^(1 I y I A) "I" ( I y I A) ~ ^n,\ tinder //j 

.-1 (2-19) 

- i5)2v/i+; + = S„j under //„ 

k=0 k=Q 

where sgn(.x) is one for jr>0 and -1 for jr<0, / is the largest integer no larger than 
(/ + y)A and 5 = (J + y)A — / Note that 1=^0 and 0 < (5 < 1 . 

Defining the per-chip signal-to-noise ratio as 

aIt, 

SNR = (2.20) 

Zi\o 

we have 




^ ^ r 4 -(SNR)-^«,i under//, 
2ol X ^(SXR)Sjo under //o 



( 2 . 21 ) 



The acquisition system does not know the exact values of S„,o and S„,, in advance; 
however, we can use nominal worst-case values for these parameters when designing the 

n—l 

thresholds. Simulation results on the schemes suggest that modeling s: 0 un- 

n-1 n-1 __ A=0 

der //, and modeling Xc*c**/and under achieve the desired decision 

A=0 A=0 

errors. Using these approximations, we have 
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under H 



( 2 . 22 ) 



under H, 



0 



The cumulative distribution function (cdQ of can be written in terms of the Q 
function as 



0 



. v„ < 0 






(2.23) 



where the Q function is defined as [Ref 8] 



QiU) 



poo 

xe~ ^ ^/q(cjc)(^jc 

3 



(2.24) 



An iterative algorithm for calculating the Q function is given in [Ref 9 ]. 

We now illustrate design techniques for determining a good choice of thresholds to 
use in the FSS test, the SPRT and the TSPRT such that our decision errors are less than 
our desired limits. We denote a and p as the desired false-alarm and miss probabilities 
respectively. 

First consider the FSS test. Using the Q function, we can write the false alarm and 
miss probabilities as 

Pfa — Qi\l r\j ) (2.25) 

P miss ~ * ~ Qi\/ »\/ ) (2-26) 

From these equations, the values of t' and M can be obtained by iteratively solving the 
equations simultaneously, such that < a and < p. 

For the SPRT we can use (2.14) to calculate the thresholds by inserting equalities 
in place of the inequalities. 

To devise thresholds for the TSPRT, we split this test into two parts: a SPRT with 
thresholds A and B and errors a,,„ and , and a FSS test with sample size n, threshold 
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T and errors and /?/„. It can be shown [Ref. 10] that the errors of the TSPRT are 
bounded bv 



^tsprl — ^spri "b ^/ss 




(2.27) 



So, to ensure < a and < /?, we split these errors into the sums of the errors 
due to the SPRT and the errors due to the FSS test. Specifically, we let 



where and p, are constants in [0,1]. By setting = p^o. and = p^P, we can design 

A A 

the thresholds A and B according to (2.14). Similarly, by setting a^„ = (l— /?„)« and 
/?/„ = (1 — /’i)/?, we can design t and /? using (2.25) and (2.26). Note that if we set 
A = A = 0, the TSPRT becomes the FSS test. If we set p^ = P\ = 1, the TSPRT becomes 
the SPRT. For values of p^ and p^ in (0,1) the TSPRT can be considered a mixture of 
an SPRT and an FSS test. 

D. PERFORMANCE COMPARISONS 

We now present simulation results showing the expected sample size and power 
function for each of our three decision schemes. The expected sample size, called the 
average sample number (ASN), is the average number of chips that the test runs before 
deciding between //„ and //,. (Note that for the FSS test the ASN is predetermined and 
equal to M). The power function is the probability of accepting //,. For all simulations 
in this thesis, the updating of the local m-sequence generator is half a chip (i.e., 
A = Vz), the value of y used is 0.5, the m-sequence used has a length of 1023 with prim- 
itive polynomial 1 -1- -f- -h w* -h w‘®, and the values of p^ and (the mixture weights 

for the TSPRT) are both 0.5. 

Figure 4 on page 18 depicts the ASN for the FSS test, SPRT and TSPRT with 
a = fJ = 0.01. Figure 5 on page 18 depicts the power functions for the same three tests. 
All functions are plotted in terms of \ J + y |. When |y + y | = 0.5 we have hypothesis 
//, (since j = 0 and y = 0.5). When \J -b y | = 2.0 we have hypothesis . The sample 
size M of the FSS test is 258. For all tests the resulting errors are equal to or slightly 
less than our design values chosen for a and /?, so our threshold modeling (2.22) works 



^ ^-sprt "b ^fss Po^- "b (I Po)^ 



P = Pspn + P/ss = P\P + (1 -P\)P 



(2.28) 
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Figure 4. ASN for Acquisition Schemes (p0 = pi = 0.5) 




Figure 5. Power Functions (p0 = pi = 0.5) 
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rather well. Under and H^ the SPRT has the smallest ASX. Both the SPRT and the 
TSPRT perform much better than the FSS test under both and Hy It has been found 
that the ASX of the SPRT is larger than that for the TSPRT for some range of |j + y | 
between 0.5 and 2.0. (It was found that the ASX of the SPRT can even be larger than 
the sample size of the FSS test when the errors a and P are chosen to be verv' small). 

Among all uncertainty phases to be tested in an m-sequence period, the synchroni- 
zation condition occurs only once, the condition 0.5 < |y -f y | < 2.0 occurs at one or a 
few uncertainty phases (depending on the value of A), and the non-synchronization 
condition (Hg) occurs at the remaining phases. So, to reduce the time taken to acquire 
the incoming m-sequence, we should minimize the ASX under Hg. One way to minimize 
the ASX under Hg is to utilize the SPRT acquisition scheme, however this results in high 
ASXs for the cases 0.5 < |y + y I < 2.0. It appears that the fastest acquisition will occur 
if we use the TSPRT with a large value of p, (2.28). This makes the ASX of the TSPRT 
approximately equal to that of the SPRT under Hg and results in smaller values of ASX 
under the cases 0.5 < |y -I- y | < 2.0 . Figure 6 on page 20 shows this comparison for the 
ASX function. Figure 7 on page 20 shows this comparison for the power function. 
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Figure 7. Pov>er Functions (p0= 0.3, pi = 0.9) 
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III. PERFORMANCE IN THE FADING CHANNEL 



A. GENERAL DESCRIPTION OF FADING 

In a communication system, the received signal may suffer time varying power 
fluctuations. Since these power fluctuations most often represent an attenuation in the 
received signal (as opposed to an amplification), the phenomenon is termed fading. The 
power instability can be considered a result of instantaneous amplification or atten- 
uation of the transmitted signal's amplitude. 

In our previous discussion, the parameters and characteristics of the communi- 
cations channel were considered to be fully known and time-invariant. In reality, espe- 
cially in cases where part of the communications path consists of an unguided medium, 
the channel characteristics vaiy randomly with time. Utilizing an unguided medium 
(e.g., the atmosphere or ionosphere) as part of the channel unavoidably subjects the 
communications system to the irregular variations which often occur in nature. Channel 
parameters can be altered due to rain and humidity (which absorb microwa\e energy), 
atmospheric density variations (which refract and reflect electromagnetic waves), 
mountains and other physical obstacles (which reflect or dissipate microwaves), changes 
in the ionosphere's electron concentration distribution (which affects refraction of High 
Frequency waves), as well as other factors. Semiperiodic and totally random variations 
in the propagation attributes of the channel can be viewed as altering the channel's 
transfer function. 

A Although fading arises from a variety of causes, fading phenomena can often be 
I modeled as causing multipath distortion. That is, many fading conditions can be mod- 
eled as causing several alternate transmission routes to arise between the transmitter and 
receiver. In addition to the transmission paths which the engineer incorporated into his 
design, additional unplanned transmission paths are excited randomly with time. A sig- 
nal component arriving at the receiver via one path may be out of phase with the signal 
component arriving via a second (and additional) path(s). The arriving signal compo- 
nents interfere (constructively and destructively), resulting in fading. 

I Fading channels are classified as flat fading channels and bandwidth-selective chan- 
*1 nels. In a flat fading channel all of the frequencies present in the transmitted signal fade 
^ in exactly the same manner and the communication signal can be regarded as undis- 
torted in relative shape. The received signal in such cases can be viewed as multiplied 



21 



by a random variable which accounts for the fading. One example of this type of fading 
occurs in space links passing through a turbulent atmosphere [Ref. 11 : pp. 131-132]. 
Flat fading is further categorized as slow fading and fast fading. In slow fading the 
channel variations (which giv’e rise to the signal power fluctuations) are slower than the 
longest period component in the waveform. In such cases we consider the received sig- 
nal to be multiplied by a constant (although random) amplitude coefficient. In fast 
fading, the channel variations are comparable to the period components in the 
waveform. In such cases we consider the received signal to be multiplied by a fluctuating 
‘ time function which effectively modulates the received signal. 

In bandwidth-selective channels some of the transmitted signal's frequency compo- 
nents are affected to a greater degree than others. Such channel effects are equivalent 
to the insertion of a bandpass filter into the signal transmission path [Ref 11 ; p.l32] 
which causes the signal's edge frequencies to have a larger or smaller fade than the center 
frequencies. In such cases a sudden attenuation of just the high frequency components 
may occur, or, alternately, a loss of only low frequency power may be observed. 

It should be noted at this point that fading is a major communications obstacle 
which can cause signal attenuations of several tens of decibels. Figure 8 on page 23 
shows the median duration of fast fades as a function of fade depth for a 4 GHz signal 
over a path length of 30-35 miles. 

In the remainder of this paper we narrow our focus to consider on\\ flat, slow fading 
channels. 

^ Many methods have been devised to counter the effects of fading. The performance 
reduction that characterizes fading is caused by the received signal's amplitude being 
weakened in comparison to its design value. One technique for minimizing the received 
signal's amplitude fluctuations is to incorporate automatic gain control (AGC) into the 
receiving amplifier [Ref 12 : p.ll2j. This is a negative feedback method from control 
theor>' in which the receiver amplifier output signal is employed to adjust the amplifier's 
gain. The circuit is designed so that the output level remains steady notwithstanding 
changes in the input amplitude. The major disadvantage of this technique is that the 
channel additive noise amplitude is also increased along with the signal amplitude. 

A better technique to counter fading consists of dividing the communication's signal 
power between multiple subchannels which fade independently of each other. The 
probability that fading will be extreme in all subchannels simultaneously is very low, so 
if all subchannel outputs are used to reconstruct the signal, we expect better perform- 
ance than over a single channel. The use of multiple transmission subchannels is called 
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Source: Pierce, J.R., and Pasner, E.C., Iniroduciion To Communication Science And 
Systems, Plenum Press, 1980 



diversity transmission (Ref. 13 : pp. 346-350]. Diversity transmission methods fall under 
three general areas (Ref 14 : pp. 632-634): 

1. space diversity utilizes several receivers in separate geographical areas. The receiver 
outputs are collected to reconstitute the desired signal. When the signal at one 
receiver fades it is hoped that the signal at a distant location (i.e., at the other re- 
ceivers) will be unaffected. 

2. time diversity employs retransmission of the same signal at spaced time intervals. 
It is hoped that if time-varying fading affects a signal at one time, it will not affect 
the signal at a later time. 

3. frequency diversity employs multiple frequency channels to combat the effects of 
frequency selective fading. It is hoped that if the signal suffers fading in one fre- 
quency band, other frequency bands will still be usable. 

B. MATHEMATICAL DESCRIPTION OF FADING 

Thus far we have described fading phenomena in general qualitative terms. Al- 
though we have asserted that in slow, flat fading channels the received signal is effec- 
tively multiplied by a constant (albeit random) amplitude coefficient, we have yet to 
grant mathematical justification to this statement. In order to proceed any further with 
the analysis of our receiver in a fading channel, we must determine a suitable math- 
ematical model to describe fading. Since fading phenomena largely depend on unpre- 
dictable variations in the environment, we can not model fading deterministically. 
(Indeed, if it could be modeled deterministically we would then simply account for it by 
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adding an amplification factor in our receiver design and no further consideration would 
be necessar> .) The best that can be achieved is a statistical description of the amplitude 
coefficient. 

Communication schemes are typically constructed to achieve a preset minimum 
standard of performance. For example, we may design a system so that the probability 
of false receiver synchronization is 0.01 and the probability of failing to detect synchro- 
nization is 0.01. Using a statistical description of fading we can simulate our system in 
a fading channel and determine if our minimum standards of performance are attainable. 
If the random channel attenuation factor degrades our system to the point where we fail 
to satisfy our performance standards we can employ three corrective options. 

Our first option is to attempt to modify the receiver design by suitably changing the 
decision thresholds to meet the performance standards. 

Alternately, a second course of action would be to increase our transmitter power 
to a level where it compensates for a worst-case power loss in the fading channel. In this 
option we meet our performance standards by allowing adequate "design margins." 

Our third approach would be to relax our performance standards. This option may 
often be the most attractive if we consider other important design variables. As an il- 
lustration of how the third approach may be the most advantageous, consider again the 
system where we desire that the probability of false receiver synchronization and prob- 
ability of failing to detect synchronization both be 0.01. Suppose further that a simu- 
lation shows that fading reduces both these specifications to 0.1. In our fading channel, 
with our degraded performance standards, it may be the case that the expected time to 
acquire synchronization is one minute. If we employ option one above and redesign the 
receiver, it may be the case that the performance standards are met, but the time to ac- 
quire synchronization is increased to one hour. In such a scenario we would probably 
sacrifice our performance standards for the sake of acquisition time. 

Consider a multipath channel. In addition to the transmission path which the en- 
gineer incorporated into his design, additional unplanned transmission paths of various 
electrical lengths are excited. Portions of the power transmitted at a given time arrive 
"sequentially" at the receiver over each of the sundry transmission paths. Since we are 
restricting our consideration only to flat fading channels, it suffices to examine the 
channel effects on a single frequency component; all other frequencies will be affected 
in a similar way. 

A particular frequency component is transmitted at a certain time. The received 
signals will arrive at the receiver over the various paths. Since the lengths of the paths 
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difTer, the received sine waves will (presumably) be out of phase with each other, and 
when combined they will constructively or destructively interfere. The fading channel 
causes the electrical path lengths to var>’ continually and randomly and this in turn 
causes the phase differences between the component tones to vary. The superposition 
of all the tones (i.e., the received signal) will have an envelope and phase that will vary’ 
randomly as a result. So, in a fading channel, the received signal is the superposition 
of a number of random phasors. 

Initially assume that each of the random phasors is of a comparable magnitude. 
Consider each of the individual phasors to be resolved into quadrature components. 
Since the fading channel randomly varies the path lengths, the quadrature components 
of one particular random phasor (at a particular time) will be independent of the 
quadrature components of all of the other random phasors. Additionally, the 
quadrature components of any particular phasor are uncorrelated (since the phasor is 
resolved into sine and cosine components). 

Before finding the total received signal, first find the total signal in each of the 
quadrature directions. We do this by simply adding the particular quadrature compo- 
nents of all the random phasors. In doing this we are adding a number of independent 
random variables (from the same underlying distribution) together, and the sum, by the 
central limit theorem (Ref 15 ; p.210], approaches a Gaussian random variable as the 
number of contributing phasors grows large. Thus, at any given time, the two 
quadrature resultants are independent Gaussian random variables. The envelope of the 
received signal will therefore be Rayleigh distributed [Ref 16 : pp.101-103). (Note that 
the phase of the received signal is of no concern since we utilize a noncoherent detection 
scheme in the receiver.) 

The received signal, considering all frequencies, will thus be effectively multiplied by 
a coefficient, ij/, having the probability density function (pdf) 

= (3.1) 

a 

where determined by the nature of the channel, is the variance of each of the 
Gaussian quadrature random variates. 

Since the above discussion utilized the central limit theorem, the question naturally 
arises: how many alternate transmission paths are required for the Rayleigh model to 
be valid? Surprisingly, it has been noted [Ref 17 : pp. 348-351] that as few as six sine 
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waves with mdependently van'ing random phases will combine to give a fluctuating re- 
sultant whose envelope closely approaches Rayleigh statistics. See Figure 9. 




Percentage of time that signal level exceeds abscissa value 

Figure 9. Envelope Distributions For n Equal-Amplitude Random-Phase 
Phasors: Solid Lines are Rayleigh Distributions. 

Source: Schwartz, M., Bennett, W.R., Stein, S., Communication Systems and Techniques, 
McGraw Hill, 1966, p349. 



In deriving the Rayleigh channel results, we assumed that each of the random 
phasors had a comparable magnitude. It is more realistic to assume that one of the 
paths, namely the path that we intend to design for (i.e., the no-fading path), is much 
stronger than the other paths. The power in the strong, intended path is termed the di- 
rect component. The total power in all the weak paths is termed the diffuse component. 
The net received signal v,ill then be a combination of a steady tone (again, considering 
a single frequency since the channel is flat) and Gaussian quadrature random variables. 
The envelope statistics for "specular-plus-Rayleigh" fading, called Ricean fading, are 
given by [Ref 17 : p.372]: 



M'l') = 



± 

2 

a 







(3.2) 



where 
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= instantaneous value of the fading coefficient 
= power in the direct component 
2a^ = power in the diffuse component 

and with the total power having a normalized value of 1. 

c2 

Letting r = we find that 

MiP) = 2^f^(l + <^^‘+^)4(2^7r(I + r)) (3.3) 

Note from (3.2) that when the ratio -::p^ approaches zero, the pdf approaches the 
Rayleigh pdf (3.1). \^Tien /• -+ oo , we have no fading. Figure 10 shows the Ricean 
probability density function (3.3) for various values of r. Note that for large values of 
r the distribution of the fading coefficient is concentrated around unity. 
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C. TEST STATISTIC DENSITY FUNCTION 

We now determine the probability density function of the receiver's test statistic in 
the presence of a slow, flat fading channel. Since we assume slow fading, ip does not 
var>‘ much within one data bit interval, and thus can be assumed constant during the 
acquisition process. 

The received signal through the fading channel is given by 

f r{i) = Ao^a{i + iATc)cos{(o^t + 9) + n{i) (3.4) 

where i/> is the fading random variable with a Ricean probability density function given 
by (3.2). Referring to our receiver structure (Figure 1), equations (2.2) and (2.3) become 

A',,„ = ^il^T,S,cos9 + .\\„ (3.5) 



i{/T,S„ sin 9 + (3.6) 

where and are still defined by (2.4) and (2.5). As before, the test statistic for 
determining alignment is 

Y„ = xl + Xl, (3.7) 



The conditional probability density function ofy»„ given ip is noncentral Chi-squared: 



/Cv’J^) = 



,-6, 






2at 



)(\/ ^■n'l'^yn I }'n ^ 0 



(3.8) 



where, as with no fading, we have — (/I^'4)7]S^ and al = nT^NJA. 

We have readily determined the test statistic's conditional density function given 
We are interested in the density function without conditionality'. To determine 
f(y„), we "integrate out" the ip dependance: 



Rvn) = >A)/(>A¥>A 
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= "'"t \J-^)diP (3.9) 

2a„a Jq \ a„ J \ a / 



To simplify notation, we combine terms by letting 



a = 



5 , I 

’ ^-~2 ’ + 

a lo^ lo 



(3.10) 



and 



d = - 



1 -s^2e^ 

e ^ " e 



> 2 2 
la„a 



(3.11) 



Substituting (3.10) and (3.11) into (3.9), we can rewrite the probability density function 
as 



RVn) 







(3.12) 



The Bessel function is defined as a power series expansion. Substituting the definition 
of the Bessel function into (3.12) yields 



/»oo 



fiyn) = d 



4/e 






Vl! 



•^0 



^ 4^ 0- 
2=0 



oo 

4/^^a^^ V 

4^0!)^ L 



A=0 



A\k\f 



d4> 



Since integration and summation are both linear operations, the above formula can be 
rewritten as 



fb'n) 



OO oo 




J=0 k=0 



_b^ 

4\k\f 



poo 

^2(/>A)+l 






dip 



It has been derived [Ref. 18 : p. 271] that 



(3.13) 
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(3.14) 



2n+\ -ax , 

X e ax 



n\ 



Applying (3.14) to (3.13) results in 



/C>'„) = d 




2J,2k 
a b 



{j+k)\ 



j=Q fe=0 






which is rewritten as 



fiy-n) = 





oo 



A :=0 



b'^‘'{j + k)\ 
4*(A:!)‘c* 



fiy'n) = 




(3.15) 



The term (j+kj! = (j+ l)(j+ 2)(j+ 3)...(j+ k) can be written using Appell's symbol 
[Ref. 19 : p. 7] as {J+ I,k). Similarly, using this notation, we can rewrite k! as (I.k)- 
Adopting this notation, the density function in (3.15) is written as 



/Cv’J 



^ OO 

d V( a^y \ V {j+l,k) [b^jAcf 
2c /_j\ 4c ) j\ (1,^) 
y=0 k=Q 



(3.16) 



The second summation is in the form of the confluent hypergeometric function [Ref. 20 
: pp. 298-320): 





(/•+ 1.-^) 



|F,(/ -1- 1 ; 1 ; b^jAc) 



(3.17) 
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Substituting (3.17) into (3.16) yields 



y=0 

To further simplify, we express jF, in terms of a Laguerre Polynomial [Ref. 21 : p. 430] 
as 



,f, (c + rt ; c ; jr) = *( - jr) 

Substituting (3.19) into (3.18), we find 



(3.19) 



y=o 

where L“(.) = Lj{.) is the Laguerre Polynomial of degree j. Using an identity involving 
Laguerre Polynomials (Ref. 21 : p. 316]: 



OO 




A:=0 



we can rewrite (3.20) as 



/W - (3,21) 



where Jq{.) is the Bessel function of the first kind of order zero. Using /= J— 1 we have 

/W= (3.22) 

But Ja{ix) = /o(x), and so we rewrite the probability density function/OJ as 
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(3.23) 



/(>„) = 




Substituting the values of a. b, c and d from (3.10) and (3.1 1) back into (3.23) we have 
our final formula for /(>■„)• 



Rv.) = 



l + r 



2cr^[(l +r) + A„l2al'] 



exp 



- [(1 + + r(;./2a^)] 

( 1 + r) + ^„'2a^ 



2jril+r)(^’„!2al){;j2al) 

1 + r + X„!2al 



(3.24) 



This interesting result states that the probability density function ofy; is still noncentral 
Chi-squared. The density function is the same form as (2.8) but the parameters and 
al are replaced by new values; 



■" iq- r 



1 + r + ?.„‘2o„ 
l+r 



The cumulative distribution function corresponding to (3.24) is 

CVJ = 1 - (? 






1 + r + ^.„'2a„ 



(1 + rXv„'a„) 

1 + r + X„,'2ol 



(3.25a) 

(3.25^) 



(3.26) 



D. SIMULATION TECHNIQUES 

Now that we have a precise description of the receiver test statistic's density func- 
tion, we wish to evaluate the receiver's performance in a fading channel. Three receiver 
characteristics are of primar>’ interest for each acquisition scheme: the actual probability 
of false alarm (PFA), the actual probability of detection (PD) and the number of chips 
needed to make a decision (called the average sample number (ASN) ).2 A fourth char- 
acteristic, acquisition time, which is the average amount of time needed to correctly 



2 Note that the actual probability of false alarm (PFA) and actual probability of detection 
(PD) are distinguished from the desired probability of false alarm (a) and desired probability of 
detection (1 — /?). 
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phase synchronize the incoming and locally generated m-sequences, depends on the 
PFA, PD and ASN. 

There are two general approaches that one may use to determine the PFA, PD and 
ASN for a given fading channel. The first approach is mathematical analysis. This 
method entails determining an equation or formula which allows one to directly compute 
the quantity of interest in terms of other known quantities. Needless to say, closed form 
analytical solutions for the receiver under consideration have eluded, and continue to 
elude, researchers. (The problem, in fact, has not even been solved for the mathemat- 
ically simpler case where the correlators are reset to achieve independent identically dis- 
tributed samples.) 

The second approach used to determine the PFA, PD and ASN for given fading 
conditions is simulation. With this method, the threshold values for the various tests are 
stored in computer memory', and a computer program is used to simulate the receiver 
operation and the received signal. When the simulation shows that a threshold is ex- 
ceeded, we stop the run and measure the ASN. By running the tests a large number of 
times under various conditions, we can confidently determine the PFA and PD. All re- 
sults presented in this thesis for ASN are based on simulation. 

To illustrate the simulation techniques, we will focus on attempting to determine the 
average ASN (for a particular scheme) with fading. The techniques for determining the 
average PFA and average PD are identical. 

We are interested in determining the effect that fading has on the detection scheme's 
ASN. The average sample number depends on many factors (e.g., noise power spectral 
density, SNR, etc). To isolate the fading effect, we hold all factors stationary except for 
the fading magnitude. So, in other words, we are interested in determining the effects 
of fading on ASN w'ith all other factors remaining constant. 

With the above restriction, the ASN is a function of the fading random variable \j/, 

ASN = gW (3.27) 

The expected value of ASN in the presence of fading is thus 

E[ASN] = E[^(iA)] = rgWnm (3.28) 

•'0 

where /(^) is the Ricean probability density function describing the fading, given by 
(3.2) or (3.3). 
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This is precisely the point where the analysis stops. If we could solve the integral 
in (3.28) we w'ould have the answer we seek. How'ever, since the ASN function (3.27) 
is presently unknown, the integral is not solvable and we must resort to simulation. The 
only purpose in performing the following simulation is to solve (3.28) for E[ASN] (and 
to solve similar equations for E[PD] and E[PFA]). 

1. Monte Carlo Simulation 

We wish to determine a value for the integral in (3.28). The strong law of large 
numbers (Ref. 22 : pp. 88-90] states that for "large enough" n, we can determine 
E[g(i/')] (and hence (3.28)) from a sample mean; 



n 

(=i 

So, our algorithm for computing /“^(•/')/(i/')^<A is very simple [Ref. 22 : pp. 181-194]; 

1. Generate n independent samples of a Ricean random variable whose density 
function describes the fading. 

2. Determine g{\p), the ASN for a particular \p„ either directly or by simulation for 
!</<«. 

1 ” 

3. The desired result is X Si'P,)- 

1=1 

Step I. above, generating independent samples of a Ricean random variable, is 
of critical importance. If the samples are not Ricean and'or not independent the algo- 
rithm will not be successful. The method used for generating such independent Ricean 
random variables is included in the Appendix at the end of this thesis. 

We characterize the severity of fading by specifying a value of r, where r is the 
ratio of the power in the direct component to the power in the diffuse component and 
the Ricean density function is described by (3.3). For each Ricean random variate ijj, 
we run a simulation program to determine the ASN of a test given that value of if/,. 
We run the simulation program for a large number of independent variates \p„ and 
compute the average ASN by the summation in step 3. above. 

The Monte Carlo method gives us only an approximation of the integral we are 
trying to evaluate. The sample mean that we compute in the third step of the algorithm 
is an estimate of the integral which becomes exact only as n oo. We are interested in 
determining the mean square error of our estimator. 

Recall that our goal is to evaluate E[ASN] given by 
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rgwfwdii^ = 



(3.30) 



where we designate the exact (albeit unknown) value of this integral by I. Designate the 
estimate of the integral by g„. That is, 






(3.31) 



1=1 



Then the mean of the square of the error of our estimator is given by [Ref 22 : pp. 
194-196] 



E[(sample mean error)^] = E[(^„ — 1)^] = Var(g„) 



(3.32) 



But by using the definition of g„, the variance of g„ can be rewritten as 
Var(gJ = Var|-^|jg(.A,)| = Var||j g(i/.,)| 



(3.33) 



Substituting (3.33) into (3.32) and using the shorthand notation 
V = E[(sample mean error)^], we have 



V = -YVar(^(iAi) + si^T) + ••• + si'kn)) 



(3.34) 



Since the variance of the sum of independent random variables is the sum of the vari- 
ances, (3.34) is revised; 



V = -V(Var[g(.A,)] + Var[g(t/.2)] + - + Vark(iA.)]} 



(3.35) 



But each of the variance terms are equal to each other. That is to say, VarQg(i/',)] = 
VariXiAz)] Var[^(iAJ] = Var[g(t^)] . Using this fact 



V = — {rtVar[^(iA)]} = — Var[g(t(^)] 



(3.36) 



But the variance ofg(i/^) is equal to E[^^(iA)] - (E[g(i/')]}^ = E[g^(t(^)] - P . Thus 
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V = - P} 



(3.37) 



But we know that 






<*00 

g\^)fWd^ 

*'o 



so we substitute this quantity into (3.37) and arrive at 



V 






n 



(3.38) 



Xow, neither of the quantities in the numerator of (3.38) is known since we do not have 
an expression for g(iA). Thus, we again resort to the law of large numbers to approxi- 
mate each of the numerator terms [Ref 22 : pp. 196-197]: 



g WfWd^ 




(3.39a) 



• ( 3 . 39 ^) 

/=i 

So, our final equation for the mean square error for our .Monte Carlo estimator is given 
by 



V = 




(3.40) 



What is so nice about the .Monte Carlo method is that as we compute our estimator we 
can simultaneously compute just how good or bad our estimator is by using (3.40). 
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Again, the exact value of the integral that we are trying to determine is desig- 
nated as /. We would like to be able to place an upper and lower bound on I with a 
certain degree of confidence. Mathematically, we would like to be able to specify a 
number d and then compute two constants c, and Cj such that 

P[c, < 1 < C 2 ] > \-d (3.41) 

Chebyshev's Inequality [Ref 15 : p. 105] states that 

P[ /z — /cr < A' < ^ -f /(7 ] > 1 ^ (3.42) 



where n is the mean of the random variable X and a is its standard deviation. Applying 
the inequality to our scenario, we have 



< I < g„ + /vV ] > 1--^ (3.43) 

Thus, to choose bounds for (3.41), simply let 

/ - -j=- (3.44a) 

Ci = gn- ^/V , C2 = + ijV (3.44^) 



It is known that Chebyshev's Inequality gives a bound that is often pessimistic 
[Ref 22 : p. 197). A much tighter bound can be obtained if we use a "large enough" 
number of random variates so that the central limit theorem is applicable. In such 
an event it can be shown [Ref 22 : p. 198] that 



Vv 



N(0,1) 



Y 



(3.45) 



where N(0,1) is the standard normal random variable which we designate as Y. If we 
choose i such that 

l-Fy = P[T>0 = y (3.46) 

then rearranging (3.43) results in 
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P( li-^J </Vv ) = p 



1 -Sn 

Vv 



< I 



( » -Jn 

V vv 





1 -d 



(3.47) 



As an example of how to apply the above method, suppose we want to deter- 
mine an upper and lower bound for 1 with 90% certainty. Using (3.41) we would choose 
d = 0.1. Using this value of d, we determine i (from tables) such that 

P[ P > O = 4 



where }' is the standard normal distribution. Now, with this value of /, we can be 90% 
certain that 1 lies between g„ — and g„ -f- r^/v" where g„ and v are determined from 
(3.31) and (3.40) respectively. 

2. Numerical Integration 

The .Monte Carlo method is but one way to evaluate E[.-\SN] given by the in- 
tegral in (3. 28). A simpler and more direct technique involves estimating the integral 
by using a numerical method. 

The integral we wish to evaluate is given by (3.30) and repeated below: 



Too 

(3.49) 

->0 

We have a closed form analytical description of given in (3.24), but we do not 
know the form of g{il/). Note that for all ij/, g{<P)f{*P) > 0. 

I Our approach is to view the integral as an area. We start by dividing the density 
function of /(lA) into a certain number of equal areas. For illustration, we will divide the 
density function into 25 sections, each with an area of 1/25 = 0.04 square units. So the 
first area spans the abscissa from 0 to The second area spans the abscissa from lAi 
to lAz, and so on. The last area spans the abscissa from iAj 4 to oo. So we can write (3.49) 
as 25 separate integrals: 

poo piAi flJ/2 f«> 

g{^)f{ii¥^ = g{^)f{^)<^'li + gWf{*P)d'P -)-•••+ g(iA)/(iAW (3.50) 

*’o •’o •'(Z'24 
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Note that the values ... easy to determine since f{\j/) is known. Also note that 

(3.50) is still exact. 

Let the midpoints of each of the 25 intervals (0,i/^,), (lAi.t^^) ••• be desig- 

nated in general as \p" and specifically as \pj^, \p^, \pc, ••• <Ay respcctively.3 Note that the 
majority of the mass of f{p) is concentrated in the region between ij/ = 0 and approxi- 
mately \j/ = 3 (depending on the value of r); sec Figure 10. Because of this, we expect 
the 25 points iAa-iAb. ••• Py to be clustered close together in the region (0,3). Putting this 
another way, we expect each of the first 24 intervals ( 0 ,i/^,),(^„i/^ 2 ) ... ( 1 ^ 23 ,^ 24 ) to be very 
"thin”. 

Now we must make an approximation. Suppose over each of the 25 areas, we 
assume g{il/) = g{\p') where is the midpoint of the interval (^,,^,* 1 ). Since ip undergoes 
a very' small variation between p, and we assume g{p) also undergoes only a very 
small variation over the interval, and is, essentially, constant with a value ofg(^’). With 
this approximation, (3.50) can be rewritten as 

roo rC '2 foo 

giP)f{P)^P-g{pA)\ fWdp+gipB)] fWdp + -+giPy)\ f{p)dp Ci.51) 
*'0 Jo J ^>24 

But, by the way that we chose p 2 ... P 24 , each integral has an area exactly equal to 
1/25. Thus we have 

\^gWfiP)dp -^IgipA) + gipB) +•••+ ^(<Ay)] (3.52) 

Our numerical integration algorithm can be summarized as 

1. Divide the probability density function f{p) into X equal areas. Determine the 
midpoint of each area. 

2. For each of the area midpoints, use a simulation program to determine g(p'). 

3. The E[ASN] is approximated by *^21g(<A’) 

A ,=i 

A major disadvantage of the numerical integration technique is that we do not 
have bounds on the estimator (or its error) as we did with the Monte Carlo method. 

3 Note that aU of the midpoints are readily computed except for py Since the 25th integral 
is improper, technically we should have py = 00 . To resolve this, we truncate the tail of the 25th 
integral when it goes below a certain very small value, perhaps 1 x 10'^, and then use the mean 
value theorem for integrals to find an "effective" value for py. 
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Most numerical integration schemes divide the abscissa into equal areas, and thus have 
well known error bounds (Ref. 23 : pp. 285-312]. In our technique the abscissa is spe- 
cifically noi divided up into equal areas, precisely so that we may solve each of the 25 
integrals in (3.51). 

Note that the numerical integration technique is not too unlike the Monte Carlo 
method. For our numerical integration scheme, the points are biased around the mean 
of the density function, as we would expect them to be if they were chosen randomly. 
More precisely, because we divided the density function into a certain number of equal 
areas, each of the ij/' in the numerical integration scheme are equally likely, and thus we 
can well imagine that they were generated by our random number generator in the 
Monte Carlo method. 

E. RECEIVER PERFORMANCE WITH FADING 

We now use the Monte Carlo and numerical integration algorithms to examine the 
effects of fading on the receiver acquisition schemes. We begin by e.xamining the per- 
formance of the TSPRT scheme in a fading channel characterized by r= 1. Note that 
this is rather severe fading; there is just as much diffuse signal as direct signal. For our 
simulation we use po = p^=0.5 (2.28), SN’R= -lOdB, a = /? = 0.01, and a sequence of 
length 1023. 

1. Expected Results 

Consider again the per-chip signal to noise ratio defined in (2.20) and repeated 

below: 



Rewriting (2.21) we have 



SXR 



SXR 



2.Vo 




under //, 
under Hq 



(3.53) 



(3.55) 



We know that fading changes the quantities and al in accordance with (3.25). 
Therefore, under both Ilg and //, it is the case that 
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^FADING 



^orioinalI 



r 



(3.56) 



1 + /• + >.„/2a^ 



In other words, fading has the effect of reducing the SNR by a factor of 
r/(l + r + . Since AJal is always greater than zero, this multiplicative factor will 

always be less than one, and thus the SNR is always effectively reduced below its original 
value, with the exact degree of reduction dependent upon the value of r.4 

Consider again our model used in determining decision thresholds, presented in 
(2.22) and repeated below: 



We first examine H^- Without fading, our design anticipates to be ap- 

proximately equal to the per-chip SNR. With fading, the actual SNR is effectively re- 
duced, so we expect our test to become overdesigned (e.g., overly conservative) under //„. 
In other words, whereas we accounted for a certain SNR in designing for a low proba- 
bility of false alarm (see (3.57)), the actual SNR is even less than what we designed for, 
thus making our correlator output appear even more "HQ-like". So, we expect fading to 
actually improve PFA (i.e., make it smaller) and lower the ASN under Hq. We expect 
our tests, under Hq, to perform better as the fading becomes worse. 

Having said this, we quickly add that we expect the effects of fading to be only 
slight under Hq. By (3.57), a reduction in effective SNR will mean that our test is over- 
designed by some fixed amount. Since the two m-sequences are not synchronized under 
Hq, we should, ideally, detect only noise after the correlation process in the receiver. 
Since we have only partial correlation in the decision process, we expect fading to in- 
fluence the test (as described above), but only to a small extent. 

The situation is much different under //,. As with Hq, our design anticipates 
XJlal to be proportional to the per-chip SNR. However, with fading, our actual effec- 
tive SN'R is less than what our test accounts for in attempting to satisfy the desired 
probability of detection. With fading, our signal appears less "H^-Uke". So, we expect 

4 Note that for r-* oo the reduction factor becomes equal to one. This is expected, since 
r = oo represents no fading. 




< SNR 



«(SNR)(l-|y|A)' 




(3.57) 
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fading to cause a deterioration in the receiver's probability of detection, PD, for all ac- 
quisition schemes. 

Although we expect the effects of fading to be slight under we expect the 
effects to be severe under //,. Examining (3.57) we conclude that under //„ the test is 
overdesigned by a fixed amount equal to the effective reduction in SNR. Under //,, 
however, the effective SNR reduction is made worse as n increases. Specifically, apply- 
ing (3.56) we see that under the SNR is effectively reduced br a constant factor of 
( a -/[1 -i- /• -h SA’i?]). Under the SNR is effectively reduced by the factor 



which becomes smaller as « increases. Our scheme under H, becomes successively more 
underdesigned as n increases. 

2. Monte Carlo Simulation Results 

A 50 run Monte Carlo TSPRT simulation was performed under fading condi- 
tions with r = 1. After the data for all runs were collected, it was arranged in order of 
increasing \{/, the fading random variable. Figure 11 on page 43 displays ASN vs. ij/ 
under H^. Figure 12 on page 43 displays PD vs i[/ under //,. Recall that \{/ = 1 amounts 
to the absence of fading. As \j/ increases beyond ^ = 1, PD remains high and ASN de- 
creases dramatically. This is expected, since \j/ > 1 corresponds to a stronger signal (see 
(3.4)) than expected. Our test is designed to meet PD = 0.99 with no fading; if we have 
"constructive fading" (i.e., if/ > \) our signal is stronger than what we have designed for 
and therefore fewer chips are needed in the correlation process to meet the Hi threshold. 

As \j/ goes below one, the ASN starts increasing rapidly while PD starts declin- 
ing. The test runs longer in attempting to satisfy PD = 0.99, but increasingly decides Hu 
by mistake. As if/ decreases, the Hi condition increasingly appears to the receiver to be 
the Ho condition. At about i/^ = 0.55 , PD drops below 0.5 and the decision scheme starts 
deciding Hq more often than //,. As ij/ drops below 0.55 the ASN drops as the TSPRT 
decides on Ho earlier and earlier. 

Figure 13 on page 44 shows E[ASN] , given by (3.31), vs. the number of runs. 
Note that E[ASN] seems to settle down as the number of runs increases. Under the 
Monte Carlo scheme, E[ASN] would give the exact answer as n -* oo. Generally, the 
more runs used, the better the estimator is expected to be. 
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Figure 11. TSPRT ASN vs ^ under (Monte Carlo) 




Figure 12. TSPRT PD vs ^ under Hy (Monte Carlo) 
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Figure 13. 




Figure 14 on page 45 displays ASN vs \p under Hg. Notice that the change in 
,ASX over the range of \]/ is only about ten chips, a small variation compared to the 
nearly 200 chip fluctuation seen under H, (Figure 11). Figure 15 on page 45 displays 
PFA vs. \j/ under Hg. Notice that PFA does not go above the design value of a = 0.01 
until i/^2=1.5. This confirms our expectation that fading would have only a small effect 
under Hg. 

How confident can we be that our Monte Carlo results are accurate? 
Figure 16 on page 46 shows the mean square error of our ASN estimate (under //,) vs. 
the number of runs. We expect that as the number of runs increases, the mean square 
error, V, given by (3.40), will decrease, and we see that this is indeed the case. Ideally, 
F ^ 0 as the number of runs -♦ oo. After 50 runs, our Monte Carlo scheme yields 
E[,ASN] = 140 with a mean square error 70. Obviously, our results after 50 runs do 
not instill as much confidence as we would have desired. 

A confidence interval calculation is instructive. Suppose we want to determine 
bounds for / with 90% confidence. Using (3.41), we see we require d = 0. 1. Using (3.48) 
we find t such that P[F> r] = 0.05 where Y is the standard normal distribution. Con- 
sulting a table yields / = 1.645. With g„ = 140 and F= 70, (3.43) gives our desired con- 
fidence interval: 
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Figure 14. TSPRT ASN vs ^ under (Monte Carlo) 




Figure 15. TSPRT PFA vs ^ under (Monte Carlo) 
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Figure 16. 




P[126<I<154] > 90% 

In other words, we can be assured that E[ASN] under //, is within the 28 chip band 
(126-154) with 90% confidence. 

3. Numerical Integration Simulation Results 

We reran the TSPRT simulation under r— 1 fading conditions by dividing the 
Ricean density into 25 equal areas. Figure 17 on page 47 shows the ASN vs for 
Note that these results, and the results to follow, closely resemble those obtained 
through Monte Carlo simulation. Figure 18 on page 47 displays PD vs. under //,. 
Figure 19 on page 48 and Figure 20 on page 48 show ASN and PFA, respectively, vs. 
^ under Figure 21 on page 49 includes data for PD vs. \}/ for points between and 
H,. 

How confident can we be that our numerical integration results are accurate? 
Recall that there is no precise way to determine confidence intervals for our integration 
scheme. To gain insight into how stable our results are, we reran the integration scheme 
but divided the density function into 50 areas instead of 25. We found that, to the 
nearest chip, the value of E[ASN] under //, for 50 areas w’as identical to the value ob- 
tained for 25 areas. Similarly, the value of E[PD] using 50 areas was found to be equal 
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Figure 17. TSPRT ASN vs under (Numerical Integration) 




Figure 18. TSPRT PD vs ^ under Hi (Numerical Integration) 
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Figure 19. TSPRT ASN vs ^ under (Numerical Integration) 




Figure 20. TSPRT PFA vs 4' under Ho (Numerical Integration) 
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Figure 21. TSPRT PD vs tl/ (Numerical Integration) 

to the value obtained with 25 areas to within three decimal places. We conclude that 
our numerical integration scheme using 25 areas is accurate and, at this point, we 
abandon the Monte Carlo method for all future simulations in this thesis. 

We summarize our results for the TSPRT with r=\. Without fading, our ex- 
pected value for the probability of detection under //, is 0.992. With fading, the expected 
value for the probability of detection dropped to 0.762. Our e.xpected value of ASN in- 
creased from 138 without fading to 156 with fading. 

Under //q our expected value of probability of false alarm decreased from the 
no-fading value of 0.01 to a value of 0.007 with fading. The ASN under decreased 
from 158 without fading to 155 with fading. 

Figure 22 on page 50 shows the TSPRT power vs. |y+y| with and without 
fading. Recall that |y -f y 1 = 0.5 corresponds to //, w^hile |y-f y | = 2.0 corresponds to 
Hf). Figure 23 on page 50 summarizes the fading results for ASN vs. ly-b y |. 

4. Simulation For Various Fading Conditions 

We obtained simulation results for the TSPRT with fading for the additional 
cases of r = 0, 10 and 20. In all cases the simulation technique consisted of numerical 
integration with the density function divided into 25 equal areas. Figure 24 on page 51 
summarizes the results for ASN vs. |y -f y | while Figure 25 on page 51 summarizes the 
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Figure 22. TSPRT Power vs |y + y | 




Figure 23, TSPRT ASN vs |y + y | 



50 






Figure 24. TSPRT ASN vs |y + y | (various r) 



Figure 25. 
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results for PD vs. |^' + y \ Note that the test performs better as r increases, as expected. 
Additionally, it is seen that fading has very little effect on the //j hypothesis, as expected. 

We then simulated the FSS test under identical conditions for r = 0,1,10 and 20. 
The results for PD vs. iy + y | are summarized in Figure 26. (Note that for the FSS test, 
the ASN is a constant independent of fading.) 



I J -t- GAMMA I 

Figure 26. FSS PD vs Jy + y 1 (various r) 

F. THRESHOLD ADJUSTMENT 

We have seen that fading causes no deterioration of the receiver performance under 
//(!, but generates a considerable degradation under Under worst possible fading, the 
Rayleigh case, the probability of detection is reduced from 0.992 to 0.724 for the TSPRT 
and from 0.990 to 0.721 for the FSS test. 

Since we know the exact form of the test statistic's probability density function with 
fading (given by (3.24)), w^e can redesign our decision processor to account for the fad- 
ing. Designating as X„_i under //„ and similarly designating X.„ as under w^e WTite 
the likelihood ratio as 
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Regrouping terms, we can isolate all terms with the quantity yJ2a\ into a test statistic 
and compare this test statistic to two thresholds; 
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where, as before, 2„J2al and ?.„^o!2al are given by (3.57) and A and B are determined by 
the procedure detailed in Section C of Chapter II. For the FSS test A = B and our two 
thresholds are equal. 

Note that the new test statistic, as well as the thresholds, depends critically on the 
value of r. Adjusting the thresholds to counter the fading effect requires that we some- 
how measure quantitatively the value of r prior to profiting from the updated design. 
In other words, if we suspect that fading has caused our receiver performance to deteri- 
orate, we must determine the degree of fading, r, prior to adjusting the thresholds. A 
detailed study of techniques that may be used to determine the fade depth is beyond the 
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scope of tliis thesis, but the following methods, individually or in combination, might 
be considered; 

1. If the fading is due to multipath, we can employ a highly directive antenna to scan 
the received signal strength in azimuth and elevation. Using this information, we 
determine how much of the received signal power is arriving from directions outside 
the design path. 

2. If the fading is due to ionospheric effects, we can use sounding techniques to 
measure the layer alterations, or we may avail ourselves of the historical informa- 
tion already gathered concerning the ionosphere. Figure 27 on page 55 [Ref 24 ] 
shows the fade depth in dB for a typical HF channel at mid-latitudes, over a 200km 
path. We can use (3.56) to convert this information into an estimate for the value 
of r. 

3. Prior to initiating communications, we can conduct an extensive channel survey to 
measure and assemble information on signal performance under var>'ing conditions 
(e.g., rain, magnetic storms, sun spot activity, channel traffic, etc). If similar con- 
ditions arise again once the channel is in use, we will have a baseline estimate of 
the amount of fading to expect. 

Howsoever determined, we assume for the remainder of this discussion that we are able 
to gauge the fade depth and calculate r. 

Figure 28 on page 56 displays the degree to which the thresholds in (3.59) var>’ with 
fading for the TSPRT with SNR = 0.1, length = 1023 and = p, =0.5. Note that the 
"original" thresholds correspond to r -> cxo, or no fading. As r decreases, fading becomes 
worse and the test runs for a longer time prior to truncation. This makes intuitive sense 
when one considers that fading effectively reduces the received signal's SNR, and so we 
must integrate more m-sequence chips to accumulate enough signal to make a decision. 
Note that if the fading becomes too severe, the TSPRT truncation point becomes greater 
than the sequence length. 

Figure 29 on page 57 shows three tj'pical test statistic trajectories and their inter- 
action with the thresholds for the case of r = 10 . Note that for the sample path with 
4/ = 0.46 (severe fading) the decision processor runs all the way to test truncation prior 
to resolving between the two hypotheses. Note that as the fading worsens from r = oo, 
100, 20 to 10 the test truncation increase from 306 to 332, 461 and 710 chips respectively. 
The truncation point for r = 1 (not shown) is an extraordinary' length of 15,221 chips. 

We reran the TSPRT with simulated fading, but with the thresholds adjusted in ac- 
cordance with (3.59), for the cases of r=10 and r=20 . As before, SNR = 0.1, 
(x = p =0.01, po and p, = 0.5 and the sequence length is 1023. The resultant power vs. 
j -t- y, the phase offset, is displayed in Figure 30 on page 58. The power curves with no 
adjustments to the thresholds, presented in Figure 25, are repeated for purposes of 
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Fade depth in dB 



Figure 27. Fade Depth For Typical HF Link 

Source; Maslin.X., HF Communications: A Systems Approach, Plenum Press, 1987 

comparison. We see that the TSPRT is fully restored; adjusting the thresholds raises the 
probability of detection to greater than 0.99. Threshold adjustment also retains the 
probability of false alarm below 0.01, as desired. 

While we may take satisfaction in restoring the TSPRT power, there is a severe price 
to be paid in employing our decision processor. Figure 31 on page 59 displays the 
TSPRT average sample number (the number of chips the decision circuit must correlate 
prior to reaching a decision) vs |y + y |. Note that the ASN has increased drastically. 

Recall that among all uncertainty phases to be tested in an m-sequence period, the 
synchronization condition occurs only once, the condition 0.5 < IJ + y | < 2.0 occurs at 
one or a few uncertainty phases (depending on the value of A), and the non- 
synchronization condition {H^ occurs at all the remaining phases. For our TSPRT 
simulation, the condition 0.5 < |y + y | < 2.0 occurs at three uncertainty phases while the 
//fl condition occurs at 2042 uncertainty phases. To reduce the amount of time taken to 
acquire the incoming m-sequence it is critical that the ASN be minimized under With 
no adjustments to the thresholds, E[ASN] under is 158 chips. When the thresholds 
are adjusted, E[ASN] increases to 215 chips for r = 20 and 308 chips for r= 10. In 
summary, the price to pay for meeting the probability of detection is a large increase in 
E[ASN] under Similar results were observed for the FSS test: the test length in- 
creased, favoring the probability of detection while sacrificing acquisition time. For 
a=P= 0.01 and SNR= -lOdB, the required FSS sample size increases from 258 (for no 
fading) to 361, 509, 6617 and 8944 when fading is characterized by r=20, 10, 1 and 0, 
respectively. 
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Figure 28. TSPRT Thresholds Adjusted For Fading 
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TEST STATISTIC 




Figure 29. Typical Trajectories For r= 10 
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Figure 30. TSPRT PD Fading Performance 
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Figure 31, TSPRT ASN Fading Performance 
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IV. PERFORMANCE WITH MODULATION 



A. GENERAL DESCRIPTION OF MODULATION 

Up to this point, we have only discussed the acquisition of a communications signal 
that does not cariy data. All the prior analysis assumed that our received signal was 

r{i) = A^aii + /ATJ cos(co^i + 0) + «(/) (4. 1 ) 

where is the signal amplitude. a{i) is the m-sequence signal waveform with phase 
iAT^ (/ taken to be an integer without loss of generality), is the chip interval, A is the 
amount that the phase of the local m-scquencc is altered during the acquisition process, 
CO, and 6 are the frequency and phase of the carrier, and n{t) is the AWGN. 

In this chapter we e.xaminc the performance of our detection schemes in the pres- 
ence of data. By far the most commonly used modulation scheme is binar>’ phase shift 
keying (BPSK). For our discussion we assume that the data rate is equal to the m- 
sequence chip rate divided by the sequence length. In other words, we assume we send 
one bit of information for each period of the m-scquence. 

In the BPSK modulation scheme, to communicate a "\" we transmit one period of 
our chosen m-sequence. To communicate a "-1" we transmit the complement of the 
m-sequence. The m-sequence can thus be modulated by the data by using a single 
EXCLUSIVE-NOR gate, with the data bit stream and m-sequence chip stream as in- 
puts. The modulated m-sequence is then multiplied by a high frequency carrier for 
transmission. Representing the data stream as d{i), our received signal is now 

r(/) = Aod{l + iATc)a{l + iAT^) cos(o)^/ -1- 0) -I- n{i) (4.2) 

By c.xamining Figure 3 and (2.6), we see that modulation may present a considerable 
problem for our receiver. Note that in the correlation process we multiply our locally 
generated m-sequence by the incoming m-sequence, and then integrate this result. If 
we transmit the data "1" during the duration of the correlation process, then (2.6) will 
apply just as before. If we transmit a "-1" during the duration of the correlation process, 
(2.6) will be preceded by a negative sign, which is subsequently removed by the squaring 
process in the decision statistic (2.7). No matter what data we send ("-1" or 'T”), the 
receiver performs exactly as it did before, so long as the data remains unchanged during 
the correlation process. 
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The problem arises when a data boundary’ is encountered, and the data changes sign 
("-1" to "I" or "I" to "-1") during the correlation process. Under the //, condition the 
correlator output builds up in magnitude until a decision is made. When a data sign 
change is encountered, the correlator output starts building in the opposite direction, 
wiping out the previous gains. As an example, suppose we start by sending a "-1". In 
the //, condition, the quantity S, , given in (2.6), will start building "up" in the negative 
direction. (The fact that S„ becomes more and more negative is no concern since we are 
interested in the square of the quantity S„.) Now, if the data bit changes to "1" during 
the correlation process, the per-chip S„ becomes positive. This positive quantity adds to 
(and thus diminishes) the previously accumulated negative running total. 

It is true that we can utilize our communications scheme in such a way that we first 
acquire the m-sequence and then, after acquisition, start sending data. Such a method 
is, in effect, utilizing a data preamble of "1111..." to assist in m-sequence acquisition. 
We are interested in seeing how our acquisition schemes perform when data is trans- 
mitted from the outset, without a "111 1 ..." preamble. Knowledge of performance with 
modulation is useful for the following reasons: 

1. If we decide to transmit a signal without data for a certain fixed period of time in 
order to assist receiver acquisition, then it is likely that in many instances the re- 
ceiver will be synchronized for a certain time before actual data transmission com- 
mences. Suppose, for example, that w'e transmit our signal without data for one 
minute in order to allow the receiver to acquire the m-sequence. Consider that 
when two m-sequences are synchronized, such a condition can be detected in about 
one second. Suppose further that it just so happens that, by chance, the receiver's 
m-sequence is immediately in phase with the incoming m-sequence at the start of 
the acquisition process, i.e., suppose that we are immediately in the //, state and 
we acquire in about one second. Then, in such a scenario, we "waste" 59 seconds; 
our receiver is "ready to receive" but no data is sent for a considerable time. If, on 
the other hand, we send data immediately, we expect that it may take longer to 
acquire but, once the receiver is synchronized, data transfer immediately com- 
mences. Now, it is true that if we send data immediately, some data sent before 
acquisition might be lost, and the time needed to acquire the signal might be 
greater. On the other hand, more data might be sent because we do not "waste" 
channel time. Knowledge of the performance with modulation may help us resolve 
this tradeoff. 

2. The spreading waveform afrj is periodic and therefore is characterized by a line 
spectrum in the frequency domain. Multiplying aflj by the cosine carrier also re- 
sults in a quantity that has a line spectrum in the frequency domain. If we are 
tr}’ing to conceal our communications from a potential adversary, we would prefer 
that our signal had a continuous spectrum instead of a more easily noticeable line 
spectrum. For this reason, in secure DSSS communications, it is not advisable to 
send preambles of all ones or all zeros. Once we start sending data, the random 
(i.e., non-periodic) nature of the data results in a continuous spectrum for our 
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communications signal. We conclude that sending data immediately affords us a 
more secure communication signal than one that begins with a "1111..." preamble. 

3. No communications scheme is perfect, and we can expect that our scheme will, at 
some time, lose synchronization after data transfer has begun. Our system may 
be operating for a period of time, after which the locally generated and received 
m-sequences fall out of synchronization and make data recover}' impossible. 
Causes for losing lock include drift in the oscillators in the transmitter and/or re- 
ceiver, equipment failure and operator error. After such an incident has occurred, 
we must go through the acquisition process all over again. During such an event 
the transmitter will (most likely) not know that the receiver has lost lock, and so 
the transmitter will not cease sending data. The receiver will have to reacquire in 
the presence of data modulation. 

B. TEST STATISTIC DENSITY FUNCTION 

We now determine the probability density function of the receiver's test statistic in 
the presence of data modulation. We make one important presupposition: we assume 
that in ever}' case the number of chips used by the decision processor is less than the 
m-sequence length. This is the same as saying that we assume that we will encounter 
at most one data boundar}' during the correlation process. 

We assume that the test statistic under is unchanged when data modulation is 
introduced. When the sequences are out of phase we have, essentially, only noise at the 
correlator outputs and so we do not expect modulation to have much of an effect. Our 
confidence in this assumption is further reinforced by the fading simulation results which 
showed that fading had virtually no effect on the hypothesis. Under with modu- 
lation, the pdf is still given by (2.8) with SNR. Similarly, the cumulative dis- 

tribution function is still given by (2.23). 

We now consider the hypothesis. When we transmit data there will be a data 
modulation boundary' every N chips, w'here N is the number of chips in one period of the 
m-sequence. We assume that the location of the data boundary is uniformly distributed 
over the m-sequence length. Specifically, if w'e start the receiver correlation process at 
t = 0, then, if / is the number of chips integrated prior to encountering the data bound- 
ar}', / is a discrete uniform random variable on {0,1,2 ... ,N — 1}. 

For each decision scheme, the correlation process will run for a certain (unknow'n) 
number of chips prior to exceeding a decision threshold. Let the number of chips in the 
integration interval for a particular test be n. In other words, our correlator runs from 
/ = 0 and makes a decision at / = nT^. The probability that the modulation boundary' 
will not appear in the integration interval is thus 1 — Now, suppose that a modu- 
lation boundar}' does appear in the integration interval, i.e., at some point in the first 



62 



n — 1 chips of the correlation process. Let the data prior to the boundary' be designated 
da and the data after the modulation boundary- be designated i/, where d^A e {1, - 1}. 
Since the data is random, the probability that a sign change occurs at the data boundary' 
is 'A :P{da¥^di)= Yi. Alternately, the probability that no sign change occurs at the 
modulation boundary is also '/z. Combining this with the fact that a modulation 
boundary will only appear in the integration interval with probability n/N, we conclude 

Prob(modulation has no effect) = ~ ~ ^ ~ ~ 2 ^ 

When modulation has no effect, the pdf of the test statistic is, as before, given by 
(2.8) with A„/2al given by n(SA/?)(l — | y | A)^ in accordance with {2.21). 

Now, we examine the one remaining case; a data boundary does occur and, at the 
boundarv’, the data changes sign. Reviewing Chapter II, we see that the only term in the 
description of the receiver decision processor which has changed is given formerly by 
(2.6), and previously modeled by (2.19). Using the same model as in Chapter II, we 
obtain a new expression for S„ in the presence of data modulation: 
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Using the fact that d^= — d^, we find 
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As before, we assume each of the summations is approximately equal to zero, so it fol- 
lows that 

d,{2i-n){\-\y\A) (4.5) 

Since (2.21) still applies, we have 
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^ data = 4- iSNR)Sl = 4- (5A’^)(2/ - n)\\ - | y 1 A) = ( (4.6) 

2<7/i change 2<7„ 



where ).„J2al is equal to w(SiVi?)(l - | y | A)^ 

Combining (4.3) and (4.6) with (2.8), we determine the closed form probability den- 
sity function for the test statistic under the hypothesis as 
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The corresponding cumulative distribution function, using the Q function defmed in 
(2.24), is given by 
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If we choose to work with a normalized random variable z, where z =yJol , then (4.7) 
is rewritten as 
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We see from (4.9) that the value of z will depend on the number of chips, n, in the inte- 
gration interval. The expected value of z for a given value of n is 

E[z] = (l--^)(2 + 2rt(SiV/?)(l-|y|A)') 

n-\ 

+ 3j^(2 + [^i^]Ws.V«)(l-|v|Ay) (4.10) 

1=0 

The quantity E[z] is plotted vs. n in Figure 32 on page 66 

C. RECEIVER PERFORMANCE WITH MODULATION 

We now use the numerical integration algorithm to examine the effects of fading on 
the receiver acquisition schemes. For all simulations we assume SNR = -lOdB and the 
m-sequence is of period 1023 chips. Additionally, for the TSPRT simulations, we as- 
sume Po = Pi = 0.5. 

1. Expected Results 

We expect that modulation will not have a severe effect on the performance 
(characterized by PEA, PD and ASN) of our receiver acquisition schemes. When a 
modulation boundary is encountered, half the time it will not represent a data sign 
change and thus will be transparent to the receiver. 

To reach a more quantitative estimate of the effects of modulation, we can ex- 
amine the mean value of the quantity ?.J2al considering its value when a data change is 
encountered (given by (4.6)) and when no data change is encountered (given by (2.22)); 

n—] 

+ ( Iv )Z!( )'(5A^^)(l - I V I A)' (4. 1 1) 

/=o 

This equation can be further simplified as 
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Comparing (4.12) to the original value of given by (2.22), we see that 

modulation has the effect of reducing the SXR by approximately 
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For large values of N, the reduction factor is small. In the worst case, n = N and the 
SNR is reduced by 1/3, or l.SdB. 

2. Numerical Integration Simulation Results: FSS Test 

We examine the effects of modulation on the FSS test assuming the original 
receiver structure. In other words, without accounting for the fact that modulation al- 
ters the test statistic's density function, we wUl determine the receiver performance. 
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We first simulate the case where the data boundary' is encountered and repres- 
ents a change in the data sign (-1 to + 1 or + 1 to -1). If no data boundarj' is encount- 
ered, or if one is encountered and it does not represent a sign change, the modulation is 
transparent to the receiver and the resulting PD and PFA are as determined in Chapter 
II. (Recall that the ASX for the FSS test is predetermined by the desired probabilities 
of false alarm and miss: a and /?.) 

Let us examine a particular case in detail. Suppose we choose a = /? = 0.01. In 
this event, the test termination (found by iteratively solving (2.25) and (2.26)) turns out 
to be 258 chips. We then run the FSS test for each possible value of the data modu- 
lation boundary’ / (i.e., we run the test for / = 0,1,2, ... 257). 

Figure 33 on page 68 displays the actual probability of detection vs. the mod- 
ulation chip boundary under the hypothesis. We see that if the modulation boundary 
is encountered early in the integration process, the test can still reasonably recover and 
yield a relatively high PD (although significantly less than 1 — /?). Similarly, if the 
modulation boundarx' is encountered ver\^ late, the integration may have built up suffi- 
ciently to result in a relatively high PD. The worst case occurs when the modulation 
boundary is midway in the integration process. In such an event, we see from (4.5) that 
S„=:0 and the decision scheme decides on //„. 

Figure 34 on page 68 shows the actual probability of false alarm vs modulation 
chip boundarv’ under the //„ hypothesis. We see that modulation does not seem to have 
a major effect on PFA. Figure 35 on page 69 shows PD vs the modulation boundary 
for some points between and 

Figure 36 on page 69 summarizes the results for the FSS test power vs \j -t- y \ 
in a scheme employing data modulation. With no data modulation, PD is 0.990. When 
data modulation is present PD drops to 0.921. The PFA remains at the same original 
value, to three decimal places, when modulation is present. 

For our given test, with our given sequence length, the modulation effect does 
not seem ver>’ severe. As a matter of fact, simulation shows that for this particular set 
of conditions (SNR, a,P,N, etc), modulation causes the same performance degradation 
as fading with r = 10. 

We now evaluate the performance of the FSS test with a variety of values of 
a and Specifically, we wish to evaluate the effects of modulation for 

0.01 <«,/?< 0 . 20 . 

To accurately evaluate the effects of modulation on the FSS test, we must first 
ascertain how the test performs without modulation for all cases of interest; in other 
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Figure 33. FSS PD vs Modulation Chip Boundarj : 




Figure 34. FSS PFA vs Modulation Chip Boundar)’: //, 
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Figure 35. FSS PD vs Modulation Chip Boundary 




Figure 36. FSS Power vs |y + y | : Modulation 
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words we must determine baseline values of PD, PFA and test termination M for all a 
and P which satisfy 0.01 < a, (I < 0.2. Figure 37 on page 71 depicts the calculated test 
termination point (assuming no modulation) vs. the design probability of false alarm 
(a) and design probability of miss (j?). The termination point is calculated using the 
methods of Chapter II. We note that the termination point seems to be slightly more 
dependent on /? than on a . 

Figure 38 on page 72 displays the actual probability of miss, 1-PD, hereafter 
designated as PM, vs. a and p in the absence of data modulation. The test performs as 
e.xpected. For a given a, PM is approximately equal to p. For a given p, PM is ap- 
proximately constant as a varies. Figure 39 on page 73 displays the results for PFA vs. 
a and P with no modulation. 

We now run the FSS test for all cases of interest (0.01 < a,p < 0.2) in the pres- 
ence of data modulation. The test termination point vs. a and p is still as shown in 
Figure 37 since , at this point, we have not adjusted the thresholds. Figure 40 on page 
74 shows the P.M as a function of the design miss and false alarm probabilities. Com- 
pare this figure to Figure 38. Figure 41 on page 75 shows the variation of PM with a 
for a fixed value of p. The results can be explained by noting that as we relax our design 
probability of false alarm, the FSS terminates sooner (see Figure 37), and thus there is 
less chance that a modulation boundarx' will be encountered in a given integration in- 
terval. So, as we increase a, modulation tends to become more transparent to the re- 
ceiver, and PM moves closer to p. Figure 42 on page 75 shows the variation of PM with 
P for a fixed a. The curves are linear (as desired), but they are offset higher than the 
ideal characteristic, which would be a line of slope 1 passing through the origin. Note 
that the degree to which the curves are offset from the ideal curve decreases as a in- 
creases. As we increase a, the test termination point decreases and modulation becomes 
more transparent to the receiver. 

Figure 43 on page 76 displays the actual probability of false alarm vs a and p. 
Note that this curve for PFA is almost identical to the results obtained when no modu- 
lation is present, shown in Figure 39. Modulation seems to have almost a negligible ef- 
fect on performance under 

To summarize our results for the FSS test, we note that modulation worsens the 
actual miss probability (Figure 38) and has a negligible effect on the actual false alarm 
probability. The degree to which the actual miss probability (with modulation) exceeds 
the desired miss probability depends on the design values chosen for a and /?. As we 
increase a and'or p the test terminates earlier (Figure 37) and therefore modulation has 
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Figure 37. Calculated FSS Termination Point 

less of an effect on receiver performance. As we increase a and/or the actual proba- 
bility of miss more closely approaches its design value. 

3. Numerical Integration Results: TSPRT 

We examine the effects of modulation on the TSPRT, assuming the original re- 
ceiver structure. Again, we designate i as the data modulation boundar>’, and we assume 
/ is a uniform discrete random variable on {0,1,2 ... N — 1}. 

When we utilized the numerical integration algorithm to simulate modulation 
for the FSS test, we essentially divided the density function of / into 1023 equal areas. 
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Figure 38. Actual FSS Probability of Miss 

Since the TSPRT simulation program takes a considerably longer amount of time to 
execute, we can not divide the density function of / into so many equal segments. In- 
stead, for all TSPRT simulations, we divide the uniform pdf into X equal areas, such that 
approximately twenty modulation boundaries may appear at equally spaced points in 
the integration interv'al. An example will illustrate our approach. For the TSPRT with 
a = jS = 0.01 , the test is truncated after 306 chips have been integrated (if the test has 
not already ended). We divide the uniform density of i into 64 equal areas of 16 chips 
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Figure 39. Actual FSS Probability of False Alarm 

each.5 We then run the TSPRT with each modulation boundary: 8, 24, 40, 56.. .296; 19 
runs in all for this case. Our final results for PD, P.M ASN and PFA will represent val- 
ues averaged over all the modulation boundaries. 



5 Since the m-sequence length is always an odd number, we either pad the sequence with an 
extra chip, or subtract a chip, before dividing the pdf of i into a certain number of equal areas, each 
with an integer number of chips. 
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Figure 40. Actual FSS Miss Probability With Modulation 

We start by examining in detail an arbitrary case of the TSPRT. We choose a 
= 0.04 and p — 0.14. In such a communications scheme, we do not want to miss the 
Hi condition more than once in seven times and we do not want to falsely declare the 
Hi condition more than once in 25 opportunities. The truncation point occurs after 162 
chips are integrated. We examine the effects of modulation with no threshold adjust- 
ment. The pdf of / is divided up such that there are 21 equally spaced potential modu- 
lation data boundaries within the integration interval. 
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Figure 41. FSS PM vs « for Fixed P With Modulation 




Figure 42. FSS PM vs P for Fixed * With Modulation 
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Figure 43. Actual FSS False Alarm Probability With Modulation 

Figure 44 on page 77 displays the average sample number vs. modulation chip 
boundary under //,. Figure 45 on page 77 shows the actual probability of detection vs. 
the modulation chip boundary under //, . We see that if the modulation boundary is 
encountered very early, the test will run longer and attempt to decide the H^ hypothesis. 
If the modulation boundary is encountered after the correlation process is well under- 
way, PD drops dramatically. The ASN also drops since the test simply decides at an 
earlier and earlier point. The PD reaches a minimum if the data modulation boundary 
is encountered at approximately the 45th chip into the receiver correlation process. If 
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Figure 44. TSPRT ASN vs Modulation Chip Boundary: Hy 




Figure 45. TSPRT PD vs Modulation Chip Boundary 
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the modulation boundar}’ is encountered after this point, PD and ASN start increasing; 
the later the modulation boundary' is encountered, the less damage done to the inte- 
gration result that has built up prior to the boundar}’. 

Figure 45 depicts the variation of PD vs. modulation chip for the TSPRT, while 
Figure 33 shows the same information for the FSS test. The minimum value of PD oc- 
curs in the FSS test at the midpoint between chip 0 and the termination chip. The 
minimum value of PD occurs in the TSPRT at a point far to the left of the midpoint 
betw'een chip 0 and the test truncation value. This difference in the two curves arises 
from the fact that, whereas the FSS always runs the whole span from chip 0 to the ter- 
mination chip, the TSPRT veiy rarely runs the whole span from chip 0 to the truncation 
chip. In actuality, for these TSPRT values of a and /?, we found that with no modulation 
the expected value of ASX for the scheme is 95 chips. So, based on this, we would in- 
tuitively guess that modulation will cause the most damage if it occurs at the midpoint 
between chip 0 and chip 95. This is, in fact, what we observe. 

Figure 46 on page 79 shows the ASX’ vs. modulation boundaiy' under and 
Figure 47 on page 79 shows the PFA vs. modulation boundaiy. Note that, as with the 
FSS test, modulation has little effect on receiver performance under the //„ hypothesis. 

We now summarize our TSPRT results for a =0.04, P =0.14. Figure 48 on 
page 80 shows the test power vs. 1^-1- y \ with and without modulation. Without mod- 
ulation, E[PD] is 0.92. With modulation, E[PD] drops to 0.89, but this is still greater 
than our desired value for the detection probability (0.86). This surprising result seems 
to indicate that out TSPRT is both overdesigned and robust. We only "request" that the 
test perform with E[PD] = 0.86. It actually performs with E[PD] = 0.92. If we then 
introduce modulation, the actual detection probability is decreased, but not enough such 
that its value is driven below our desired detection probability. This result can not be 
overemphasized: our original TSPRT scheme performs in accordance with our design 
values of a and /? even in the presence of data modulation. Figure 49 on page 80 sum- 
marizes the expected value of ASN vs. |y-l- y | with and without modulation. The effect 
of modulation on ASN is extremely small. 

The above results, while certainly encouraging, apply only to the particular case 
of a =0.04, /? =0.14. We wish to know just how robust and overdesigned the scheme 
is for other values of a and p. 

We will pick a very stringent value of a and examine the TSPRT performance 
for various values of-;?. Specifically, we let a =0.01 and let P take on values from 0.01 
to 0.12. (Equivalently, our desired probability of detection ranges from 0.88 to 0.99.) 
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Figure 46. TSPRT ASN vs Modulation Chip Boundary’: 




Figure 47. TSPRT PFA vs Modulation Chip Boundary' 
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Figure 48. TSPRT Power vs |y + y | With Modulation 




Figure 49. TSPRT ASN vs \ j + y \ With Modulation 
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Figure 50 on page 82 shows the actual detection probability (from simulation) vs. the 
desired detection probability with and without modulation. Also shown is the optimal 
curve for the actual vs. design probability of detection. Note that even with this verj' 
stringent value of a, the actual detection probability with modulation will satisfy our 
specification if our desired detection probability is less than 0.87. From Figure 50 we 
see that the degree to which the TSPRT is overdesigned increases as we relax the design 
value of detection probability. Figure 51 on page 82 displays the actual ASX vs. the 
design detection probability with and without modulation. We see that modulation has 
a negligible effect on the TSPRT ASX. 

In summar}', modulation does not greatly affect out TSPRT sequential detection 
scheme. With a very strict desired false alarm probability of 0.01 and a desired detection 
probability of 0.99, the actual detection probability in the presence of data modulation 
only drops to 0.95. If the desired detection probability is less than 0.88, the receiver will 
still perform in accordance with specifications. Xote that this was not the case for the 
FSS test. For the FSS test, our results show that given any a and p, modulation will 
drive the receiver out of specification. (The degree to which the FSS receiver is out of 
specification does, however, depend on the values of a and /?.} 

D. THRESHOLD ADJUSTMENT 

We have seen that modulation causes no deterioration of the receiver performance 
under Hq, but does cause some degradation under Hi, particularly for the FSS test. Since 
we know the exact form of the test statistic's pdf (given by (2.8) for and (4.7) or (4.9) 
for Hi) we can redesign our decision processor to account for the modulation. We 
consider the FSS test. 

Using the Q function, we can write the FSS false alarm and miss probabilities as 




+ 




«-i 




i=0 





(4.14) 
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Figure 50. TSPRT PD With Modulation 




Figure 51. TSPRT ASN With Modulation 
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We attempt to iteratively solve these equations to determine the values of t' and the 
termination point M. Unfortunately, we find that for many values of a and /? the 
equations have no joint solution. 

To investigate this dilemma, we first examine the test statistic pdf under (with 
modulation), given by (4.9), for various values of n, the number of chips in the inte- 
gration process. We have two unknowns to determine: (or z' in our normalized 

notation) and M. Equation (4.13) determines the value of z', since z' is the only un- 
known in (4.13). This value of z' is used in (4.14) to determine M. We try consecutive 
values of M in (4.14) until we find the first value of M that satisfies the equation. 

We find that sometimes no value of M exists which satisfies (4.14) once z' is deter- 
mined from (4.13). Figure 52 on page 84 shows a plot of (4.9), the test statistic pdf 
under for various values of «. Note that the value of the density function at z = 0 
increases as n increases. 

A particular example will illustrate the problem. Suppose we choose a = ^ = 0.01 
and we desire the solution of the pair of equations (4.13) and (4.14). Solving (4.13) we 
find that z'= 10.0777. Using this value of z' in (4.14), we attempt to adjust n until the 
equation is satisfied, but we find no solution. Instead of solving (4.14), let us try' to de- 
termine a graphical solution. 

Noting that (4.14) is the cumulative distribution function for the pdf under //i, we 
can rephrase (4.14) in words: determine a value of n such that the area under the prob- 
ability density function, given by (4.9) and shown in Figure 52, from z = 0 to 
z= 10.077 is less than or equal to 0.01. It turns out that the area from z = 0 to 
z = 10.077 can not be made less than 0.01. We determined the exact value of the area 
under the pdf between z = 0 and z = 10.077 for all values of n from « = 0 to « = 1000. 
The results are shown in Figure 53 on page 85. We find that the minimum area ever 
possible between z = 0 and z= 10.077 is 0.0777 square units, which occurs when 
n = 280. Returning to (4.14), we conclude that if we attempt to adjust the thresholds, 
the best miss probability we can theoretically achieve (when a =0.01) is )? = 0.0777, and 
this occurs if we run the FSS for 280 chips. 

We draw two conclusions from our specific example. 

1. If we choose a =0.01 then, no matter how we attempt to adjust the FSS scheme, 
we can not theoretically achieve a value of p less than 0.0777. If we choose a value 
of P greater than 0.0777, say 0.09 or 0.15, we should be able to adjust the FSS 
scheme such that it performs within specification when modulation is present. 

2. The value of n obtained for the smallest P, n = 280 in our case, is the "optimal" n 
we can choose if we desire a p less than 0.0777. In other words, if we desire 
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Figure 52. Density Function With Modulation: /f, 

fi = 0.05 then, if we pick a test length of 280 chips, although we will not achieve 
our desired /?, any other value of M (less than or greater than 280) will perform 
even worse. 

Note that the example we chose above is particularly strict. For many values of a 
and ^ (other than the ones we used above) the two equations (4.13) and (4.14) will have 
a joint solution which gives us a successful adjusted FSS test. Only for very small values 
of fi do the equations not have a solution. 

Figure 54 on page 86 shows the value of M, the adjusted FSS test length, for 
0.01 < a,/? < 0.20 with modulation. The optimal test length is found by solving (4.13) 
and (4.14) simultaneously for each pair (a,)5). The flat portion of the curve, for small 
values of design miss probability, represent "best" choices for M, i.e., choices for M 
which will not satisfy our desired a and /? but are, nevertheless, the best we can do. 
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Figure 53. Area Under pdf From z = 0 to z= 10.077 

Figure 55 on page 87 displays the required test length vs. desired miss probability for 
various fixed probabilities of false alarm. We expect the test to meet the specifications 
(a and /?) for all values of the design miss probability except for points on the flat 
portions of the curves in Figure 55. For the flat portions, using the indicated value of 
M for the FSS test length will give us the best miss probability that can possibly be 
achieved, although we expect /? will not be satisfied. 

Figure 56 on page 88 shows the results of actual miss probability vs a and /? with 
modulation present and thresholds suitably adjusted. Compare these results to those 
obtained without threshold adjustment, shown in Figure 40. We see that our threshold 
adjustment does indeed work where we expect it to, and does indeed fail where we pre- 
dicted it would. Prior to threshold adjustment, we see from Figure 40 that the test al- 
ways failed. Figure 57 on page 89 shows the actual miss probability vs. a for several 
fixed values of Note that if the desired miss probability is 0.1 or 0.2, the adjusted FSS 
test satisfies the performance criteria. If the desired miss probability is 0.01, we see that 
the test fails, but the curve in Figure 57 shows the best that is achievable. 

Looking at the results another way. Figure 58 on page 89 displays the actual miss 
probability vs. the design miss probability for several fixed values of a. For each value 
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of a, the adjusted FSS test performs within specification until the design miss probability 
goes below a certain value. 

Finally, Figure 59 on page 90 shows the actual false alarm probability vs. a and P 
for the FSS test with the thresholds adjusted to account for the data modulated received 
signal. Note that the results are similar to those obtained without threshold adjustment 
(Figure 43). Again, modulation does not significantly affect the test performance under 
the Ho hypothesis. 




Figure 54. Adjusted FSS Test Length 
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Figure 55. FSS Test Length vs (i : New Thresholds 
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Figure 56. Actual FSS Miss ProbabiUty: New Thresholds 
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Figure 57. Actual FSS Miss Probability vs a : New Thresholds 




Figure 58. Actual FSS Miss Probability vs /J : New Thresholds 
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Figure 59. Actual FSS False Alarms Probability: New Thresholds 
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V. CONCLUSIONS 



In this thesis, we have studied the effects of fading and modulation on noncoherent 
m-sequence acquisition schemes. 

We have considered two schemes: the fixed sample size test (FSS) and the truncated 
sequential probability ratio test (TSPRT). .Most prior studies have assumed that co- 
herent demodulation is available during the acquisition process. We have dropped this 
unrealistic assumption since, in practice, m-sequence acquisition is generally performed 
by the receiver prior to carrier recover>\ Additionally, most prior research has assumed 
that the samples used in the decision processor are independent and identically distrib- 
uted. We do not reset the integrators during the correlation process, resulting in an en- 
hanced effective signal to noise ratio. 

The effects of channel fading on the performance of the synchronization tests were 
investigated. We considered a Ricean fading channel characterized by r, where r is the 
ratio of the received signal power in the direct component to the received signal power 
in the diffuse component. The probability density function for the receiver's test statistic 
was derived for the Ricean fading channel and found to be noncentral Chi-squared. It 
was determined that fading theoretically reduces the received signal's SNR by a factor 
of r',{\ -f r + y^.„'2a;]. 

To simulate fading we considered two simulation techniques. The first was a Monte 
Carlo algorithm derived from the strong law of large numbers. It was found that the 
Monte Carlo method did not provide a very tight confidence interval for a 50-run sim- 
ulation. We then applied a numerical integration method that provided answers which 
converged within 25 runs. 

We examined the effects of channel fading on the performance of the original deci- 
sion processor. For the FSS test we observed that fading had virtually no effect on the 
false alarm probability, but drastically reduced the detection probability. With very se- 
vere fading, a detection probability of 0.99 can be decreased by more than 25%. As r 
increases, the detection probability approaches the original value without fading. For 
the TSPRT we also observed that the false alarm probability was essentially unaffected 
by the fading, while the detection probability was severely reduced. For the TSPRT, the 
ASN under also remains essentially the same, but the ASN under //, increases as the 
fading becomes more severe. 
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VVc then considered performance of the tests designed with the knowledge of the 
fading. We adjusted the sample size and thresholds of the decision schemes in accord- 
ance with the new test statistic density function. For the FSS test the performance cri- 
teria were then satisfied but the required sample size increased to a ver>’ large number 
when the fading was severe. As an example, for a = ^ = 0.01 and an SNR of -lOdB, the 
required FSS sample size increased from 25S (for no fading) to 361, 509, 6617 and 8944 
when fading is characterized by r= 20, 10, 1 and 0, respectively. Since the //j condition 
is unafrected by fading and occurs for all but a few of the uncertainty phases, adjusting 
the FSS thresholds greatly increases the acquisition time. Our results suggest that it may 
not be worthwhile to adjust the decision processor to account for fading. 

When we adjusted the thresholds for the TSPRT test, a similar result was observed. 
The updated test performed in the fading channel in accordance with specifications, but 
at the expense of a much larger ASN under //q. .-\s an example, for a = ^ = 0.01 and an 
SNR of -lOdB, the ASN under //q increased from 15S to 215 and 308 for r= 20 and 10, 
respectively. To reduce receiver acquisition time it is critical that the .ASN be minimized 
under //j. In adjusting the test thresholds, the price to pay for achieving the probability 
of detection is a large increase in ASN under //g. We conclude that it might be prefer- 
able to apply a test designed under no-fading conditions; we may simply accept any de- 
gradation in detection probability caused by the presence of fading in the channel. 

We then investigated the effects of data modulation on the performance of the syn- 
chronization tests. Again, an expression for the density function of the receiver test 
statistic was derived. To see the effects of data modulation on the performance of the 
original decision processor, we first considered tests designed under no modulation and 
evaluated their performance in the presence of data modulation. For each scheme the 
modulation effect was found to be not severe. For example, if a = /? = 0.01 and SNR 
= -lOdB, modulation causes almost the same effect on the FSS test as fading with r = 10. 
.Modulation has virtually no effect on the FSS test false alarm probability, but does 
moderately reduce the detection probability. The degradation in the detection proba- 
bility tends to lessen as we relax the design criteria (a and fi). 

We found that modulation had even less of an effect on the TSPRT scheme. For a 
stringent desired false alarm probability of 0.0 1, the scheme will meet any desired de- 
tection probability less than 0.87 even if data modulation is present. If the desired de- 
tection probability is greater than 0.87, the scheme will perform slightly worse than 
desired. For example, if the desired probability of detection for the original receiver is 
0.99, it is only reduced to 0.95 if data modulation is present. The ASN of the TSPRT 
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is essentially unchanged with data modulation. We conclude that the TSPRT is rather 
robust in the presence of data modulation and threshold adjustment may be not war- 
ranted. 

We finally considered performance of the FSS test designed with the knowledge of 
data modulation. We found that there are values of (a,j?) for which it is impossible for 
the adjusted FSS receiver to perform within specification. For these points we deter- 
mined the best sample size which resulted in the largest possible detection probability. 
For the other values of (a,^), the test can be adjusted to perform satisfactorily at the 
expense of increasing the sample size. 
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APPENDIX GENERATION OF RICEAN RANDOM VARIATES 



To generate random variates from an underlying probability distribution, one starts 
with uniformly distributed random variates (on the interval [0,1]) and uses a suitable 
transformation to arrive at the desired distribution. Uniform variates are always the 
logical starting point: the theory is ver)’ well developed and many algorithms have been 
written and thoroughly tested. 

Obviously computer-generated "random" numbers are not random; they are, in fact, 
totally predictable. Furthermore, their predictability is often a desirable property [Ref 
25 ; pp. 132-1-15]. The best we can achieve with a computer random number generator 
is a sequence of numbers which appear to be random when subjected to standard statis- 
tical tests. Generally, the best random number generators put out sequences of numbers 
which have long periods, i.e.. which take a long time to start repeating. 

Making use of abstract algebra, random number generators of the form x„ = 
/(-Vi) have been developed which have the ma.ximum possible period prior to repeating. 
The most commonly used maximum length generator is the congrueniial generator which 
generates uniform random variates u, by using the algorithm [Ref 26 ; pp. 20-26]: 

= ( oXj -I- c) ( mod m) {A.\) 



To get the period equal to m, we need to follow three rules [Ref 22 : p. 13] 

1. c must be relatively prime to m 

2. a = I (mod p) for ever)' prime p which divides m 

3. a = m (mod 4) if m is a multiple of 4 

To optimize the statistical properties of the random variates, we adhere to a fourth 
rule [Ref 27 : pp. 155-156] 

4. the multiplier a should be larger than preferably larger than m/100, but 

smaller than m — ^jm 

The simulation programs were designed to run on a VAX computer using a base of 
2 and a word size of 31. Thus theoretically, the best choice for m (to result in the longest 
sequence period) would be m = 2”. However we can not choose m= 2’’ because the 
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largest number generated by the algorithm (A.1),(A.2) is a{m — l) + c. (This occurs 
when ax,-h c = m — I . In such a case (ax, + c)(mod m) will equal m — 1 and then at,*, 
will cause a number a(m — 1) + c to be generated.) If we choose m = 2^' then, during the 
algorithm, at some point, the computer will want to work with a number larger than 2^' 
(potentially as large as <7(2’’ — 1) + c ), which will cause overflow. 

Rule 4. above suggests a possible choice for the constants a, m and c. If we choose 
a ^ then the largest number generated by the algorithm will be <7(m— l) + c 
= m”. We set this equal to the largest number obtainable in the computer: 



So, to the nearest power of two, a good choice for m is thus m = 2-° . The reason for 
choosing a power of two is to facilitate choosing the constant c which, by the third rule, 
must be relatively prime to m. Any odd number can now be chosen for c. 

Summarizing the above results, we can satisfy rules 1.-4. and have a maximum pe- 
riod for our VAX computer by choosing: 

1. m = 2“ = 104S576. This will be the number of uniform random variates gener- 
ated prior to repeating. 

2. Since we need a > JTn we must choose a > 1024. To satisfy a = 1 (mod2) and 
a = 1 (mod 4) we choose a = 1029. 

3. For c pick any odd number, subject to 4. below. 

4. The maximum value generated by the algorithm is a(2’“ — 1) -I- c s: 2’“ -I- c. Since 
this number must be kept less than 2”, we require that the number c be kept less 
than about one billion. 

For the remainder of this discussion, suppose we have written two separate uniform 
random number generators (e.g., by applying the choices outlined above but with two 
different values of c) which output independent variates U, and Uj. We now use these 
to generate other needed distributions. 

It can be shown (Ref 26 : pp. 67-68] that using the random variable U,, we can 
generate a random variate w'ith an exponential distribution of parameter one by the 
transformation 



m 



3/2 -,31 

= 2 



(A.Z) 



E 




(T.4) 
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Now, with this exponential variate E, and the independent uniform random variable 
Uj, we can generate two independent normal variates with distribution X(0,1) by using 
the procedure [Ref. 28] 

N, = ,,/2E, cos2;rU2 (A.5a) 

N '2 = ^^2E, sin 27 tU2 (^.5^) 

Suppose that instead of standard normal variates X, and X '2 we desire normal 
variates described as X(5,o^) and X'(0,a^). That is, we desire that both of our normal 
variates be independent but with variance (instead of one) and with means s and 0 
respectively. Such normal variates can be generated by using the basic properties of 
mean and variance [Ref 15 : pp. 96-103] as 

X\ = aX, + s (A. 6a) 

Xg = 0 X 2 (A. 6b) 

We have now generated two normal variates with X\ ~ S{s,a-) and Xg ~ X'(0,o^). 
To generate a Ricean distributed random variable, we need only to calculate [Ref 1 : pp. 
30-31] 

^ = \/n1 + Nb (AJ) 

Such a random variable will have a Ricean density function given by (3.2) and repeated 
below: 







(T.8) 



Recall from Chapter 111 that fading is usually characterized by specifying the value 
of f, where r is the ratio of the power in the direct component to the power in the diffuse 
component: 



r = 




with the constraint that total power has a normalized value of 1 (i.e., -b 2a^ = 1 ). 
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Wc consider an example. Suppose that we wanted to generate independent variates 
from a Ricean distribution characterized by r= 1. One way to choose and would 



be 




= 1/2 , 


= 1/4 



so that + 2o^ =1. So starting with our standard normal variates N, and Nj, we 



calculate, using (A. 6): 




Na = 


■ + Jy 




N2 

2 


and then determine \p from (A. 7). 
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